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THE VELOCITY OF PROPAGATION OF 
ELECTROMAGNETIC WAVES! 


By J. T. HENDERSON AND A. G. MUNGALL 


ABSTRACT 


Measurements of the velocity of propagation of electromagnetic waves have 
been carried out at a frequency of 400 Mc./sec., which is close to that employed 
for geodetic surveying in Canada. The results were consistent, giving a value for 
c of 299,780 km./sec. This value is probably low because the effect of disconti- 
nuities at the cavity junctions cannot be computed with certainty. 


INTRODUCTION 


In recent years there has been a renewed interest in the determination of the 
velocity of propagation of electromagnetic waves (1-8). One of the methods of 
determination involves the measurement of the resonant frequency of a 
cavity of exactly known dimensions. In the present research, a coaxial cavity, 
excited in its principal (TEM) mode at a frequency of 400 Mc./sec., was 


employed. 

The choice of a coaxial cavity and also such a low frequency, with its 
inherently low Q and attendant low accuracy of resonant frequency deter- 
mination, was made because a determination of c at a frequency close to that 
employed for geodetic surveying in this country (9) seemed to be of some 
interest. Once the frequency had been selected, the choice of cavity type was 
severely limited, since at frequencies of this order resonant cavities become 
excessively large. Indeed, the frequency of 400 Mc./sec. was substituted for 
the 200 Mc./sec. survey frequency at an early stage of the experiment because 
of constructional difficulties involved in making cavities for the lower fre- 
quency. Also, the use of a coaxial cavity required, ideally, that only the length 
of the cavity be determined precisely, the diameters of inner and outer con- 


ductors having to be constant only. 
THEORY 
The resonant frequency, fo, of a coaxial cavity resonator excited in its first 
order TEM mode is related to its dimensions and the phase velocity of pro- 
pagation of the waves by the expression 


fo = c/2L(1+1/2Q)n 


1Manuscript received January 6, 1955. ; 
Contribution from Division of Physics, National Research Council, Ottawa, Canada. Issued 


as N.R.C. No. 3603. 
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where c is the velocity of propagation of electromagnetic waves in free space, 
L is the cavity length, and n is the refractive index of the medium inside the 
cavity. Q is the quality factor of the cavity, the ratio of power stored to power 
dissipated per cycle, and is assumed to be a function of uniform, finite wall 
conductivity only. If the value of Q is lowered by other forms of losses in the 
cavity, due for example to localized resistances at the junctions of cavity walls 
and ends, then the correction factor 1+1/2Q may not completely compensate 
for all losses. 


EXPERIMENTAL 


(A) Construction 

A solid copper cavity, illustrated in Fig. 1, was built to resonate at 400 
Mc./sec. All conducting surfaces were polished to a high degree of uniformity, 
the flat surfaces being optically ground and polished, and the cylindrical sur- 
faces honed. The diameters of the inner and outer conductors were approxi- 
mately 3.81 and 8.89 cm. respectively, and the inside length about 37.5 cm. 
This choice of diameters was made to inhibit excitation of higher order modes 
of oscillation, to facilitate accurate construction, and to obtain as high a Q as 
possible. The cylindrical conductors and end plates were sufficiently thick that 
elastic distortion was small. 

Good electrical contact at the junctions of end and cylindrical surfaces was 
found to be essential, since surface currents flow back and forth across the 
junctions from the cylindrical to the flat surfaces of the end plates. Conse- 
quently, after some experimentation, copper-plated annular steel inserts were 
fitted to the end plates, and the ends of the cylindrical conductors were 
bevelled as shown in the figure. 


(B) Measurements 


The resonant frequency and Q of the cavity were both determined from 
power response curves obtained by measurement of the output power as a 
function of frequency, input power being held constant. The power measuring 
device, a crystal whose rectified output was measured by a vibrator-type, d-c. 
microvoltmeter, was calibrated and found to be linear to about 2 or 3%. 
Linearity of the detector was important since non-linearity results in distortion 
of the resonance curve with consequent error in determination of Q. The 
constancy of input power was ascertained by measurement of the standing 
wave ratio in the input line to the cavity, and measurement of the output 
power from the signal generator. For the probe penetrations used, about 
5 mm., the voltage standing wave ratio in the input line remained constant 
throughout the frequency change of the cavity resonance curve. Effect of 
probe penetration on the resonant frequency was eliminated insofar as possible 
by retraction of the probes to a point where further retraction produced no 
measurable change in resonant frequency. 

Frequency measurement was accomplished by the use of a Hewlett-Packard 
counter operating from a standard 100 kc./sec. signal which was in agreement 
with other standard time signals to a few parts in 10’. The 400 Mc./sec. signal 
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applied to the cavity was produced by a set of frequency multipliers with a 
fundamental in the range 11.104 to 11.110 Mc./sec. The frequency of the 
fundamental could be measured by the counter to about one part in 10°. 
Output power from the multipliers through the frequency spread of the 
resonance curve was constant within 1%. 

The index of refraction of the air in the cavity was determined from measure- 
ment of air temperature, pressure, and humidity, and the use of the expression 
given by Essen and Froome (6) relating 7 to these quantities. 

As mentioned in the introduction, knowledge of the absolute diameters of 
the inner and outer conductors was not required, but variations from these 
diameters had to be known accurately. The alteration of the resonance fre- 
quency of an ideal cavity due to wall perturbations may be given by the 
following expression (10): 


eco w(1-3 fez) 


where w, is the ideal resonant frequency, w is that of the perturbed cavity, and 
the integral is that of the difference of the squares of the normalized magnetic 
and electric fields over that volume of the ideal cavity removed by the wall 
perturbation. The fact that such a cavity oscillating in its principal mode 
supports electric and magnetic fields which are rotationally symmetric about 
the geometric ‘axis of the cavity, and also respectively symmetric and anti- 
symmetric about the mid-point of the length of the cavity, means that certain 
types of wall distortions produce no change in resonant frequency. 

Variations of the diameters of the two cylindrical conductors were measured 
along four diametral planes at intervals of 2} cm., to an accuracy of about 
+3, for the outer conductor and about +1y for the inner. (Dimensions are 
given in the metric system although all measurements were carried out in the 
English system of inches.) These measurements indicated that the central 
conductor was bowed by about 2u, but that the diameter was constant roughly 
to the accuracy of measurement. 

An average diameter profile of the outer conductor indicated that the cavity 
wall tended to be pushed in near one end by about 2}y and pushed out near 
the other by about 3u. These distortions were such as to tend to cancel each 
other, and the resulting perturbation of resonant frequency was less than one 
part in 10°. 

For the TEM mode of oscillation, the electric and magnetic fields are larger 
near the inner than the outer conductor. A given diameter variation of the 
inner conductor can thus be shown to cause about twice the change in resonant 
frequency that would arise from the same variation in the outer conductor. 
For example, an annular distortion in the outer conductor at the center or at 
either end, of 3.8 cm. in length and 43u in depth, would alter the resonant fre- 
quency by about 1.2 parts in 105. A similar distortion on the inner conductor 
would result in a frequency shift of about 2.5 parts in 105. 

The uncertainties in the inner diameter measurements may cause errors in 
the resonant frequency as large as 1 part in 300,000. It would not have been 
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difficult to improve the diameter measurement, but this was not done because 
of the existence of other larger errors. 

The end plates, with inserts in place, were optically ground and polished 
flat to about }y over most of their surfaces, although an edge roll of about liu 
on one plate and about 3y on the other extended inward less than a centi- 
meter from the outer insert. A check by interferometric methods of the end 
plate surfaces after the cavity had been dismantled revealed erratic distortions 
along the inserts, amounting to several fringes (for the light source used, 
1 fringe = ju), and also a more serious permanent deformation of the over-all 
surface flatness. One end plate had been distorted into a cylinder to the extent 
of about 1 fringe, and further had an increased convexity of about 1 fringe. 
The other end plate, which had fitted the cylindrical conductors slightly more 
tightly, had been asymmetrically distorted in such a way that one half of the 
formerly flat annular surface had been bent downward so that the outer edge 
of this half was about 3 fringes lower than the other. These permanent defor- 
mations observed after the cavity was dismantled probably indicate larger 
elastic deformations in the end plates when the cavity was assembled, although 
the magnitude of such deformations is difficult to estimate. 

In order that the inside length of the cavity might be determined while the 
cavity was being electrically excited, the thickness of both end plates, taken 
as the distance from the outside surface to the center of the average plane of 
the conducting annular area, was first determined to ju. The cavity was then 
assembled and placed, one end down, on a surface plate. The outside length 
was then compared, by means of a dial indicator, with interferometrically 
calibrated gauge blocks. Different sets of gauge blocks were used, and measure- 
ments agreed within +23y, the accuracy being chiefly governed by the accu- 
racy of reading of the dial indicator, which was used solely as a null instrument. 
In view of the distortion of the end plates mentioned above, measurement of 
the outside length to a greater accuracy was not required. 


(C) Determination of ‘‘c’’ 

Preliminary determinations on a coaxial cavity not possessing the copper- 
plated steel inserts on the end plates led to values of c as much as 100 km./sec. 
below the value presently accepted. The Q of this cavity was roughly 14% of 
the theoretical value of 6700. Further experiments using the inserts, and re- 
peated plating of the latter, indicated that contact resistance at the junctions 
of flat and cylindrical conducting surfaces was the cause of both the low Q and 
the low value of c. It was found possible to increase Q to a maximum of about 
3800. Change of the contact resistance, (which could be due to oxidation of the 
copper), resulted in a time variation of Q, which was sufficiently rapid that 
alterations of as much as 24% took place over a period of about two hours. 
Dismantling the cavity also affected Q, because of wear on the contacting sur- 
faces. This necessitated replating whenever the cavity was taken apart. 

Measurements of ¢ were carried out in the following manner. The cavity 
was first placed end-down on a large surface plate. The top end plate was then 
removed and placed beside the cavity. Correctly selected length bars and gauge 
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blocks were wrung together and placed in position on the surface plate, and 
the dial indicator fixed at the correct height to a flat-based stand placed beside 
the gauge blocks. The equipment was then allowed several hours to come to 
room temperature, which did not vary by more than 0.2°C. over a period of 
two or three hours. The air temperature, pressure, and humidity were then 
measured, the top end plate was replaced on the cavity, and the outside 
length of the cavity measured as previously described. 

Measurement of the frequency response of the cavity followed immediately, 
15 or 20 points being taken in order that a resonance curve might be plotted. 
Measurement of the frequency response was then repeated several times. 
After the Q dropped appreciably, or after the temperature altered more than 
0.1°C., the top only was removed from the cavity, the contacting surfaces 
replated and the above procedure repeated. 

During this routine, the pressure and relative humidity of the outside air 
remained constant within the limits of measurement. The resonant frequency 
and Q determined from each of these curves, in conjunction with the previous 
length, and air temperature, pressure, and humidity measurements were used 
to compute a value of c, and each such computation was considered as a sepa- 
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rate measurement of c. The cavity was dismantled four times so that five samples 
of air were used. The resonant frequency was considered to be the center value 
of the resonance curve, and the Q was determined from the half power response 
frequencies determined from the curve. A typical curve is shown in Fig. 2. 

The wavelength, our second independent variable, is fixed by the internal 
length, L, of the cavity. However, L is less than a half-wavelength because the 
wave penetrates the two end walls. A simple calculation shows that this effect 
increases c by 2.7 km./sec. corresponding to an increase in the wavelength of 
6.6 u. Sixteen measurements were made using five samples of air. The means 
of the five sets, weighted as the square root of the number of measurements, 
give a mean value for c of 299,780 km./sec. The individual values are given in 
Table I. 

















TABLE I 
Corrections to ¢, 
Resonant Mechanical km./sec. 
Observa- frequency, length, Sennett os Gj 
tion Mc./sec. cm. Q Q Air km./sec. 
1 399 .859 37 .4672 3470 43.3 102.9 299,780 
2 399 .854 37 .4672 3310 45.3 102.9 299,778 
* 
3 399 .899 37.4641 3703 40.6 99.6 299,779 
*x 
4 399 .881 37 .4654 3376 44.5 99.6 299,780 
5 399 .877 37 .4654 3000 50.0 99.6 299,783 
6 399 .877 37 .4654 3000 50.0 99.6 299,784 
7 399 .877 37 .4654 3100 48.4 99.6 299,781 
* 
8 399 .872 37 .4657 2710 55.3 99.3 299,786 
a 
9 399 .870 37 .4659 3310 45.3 99.3 299,777 
10 399 .870 37.4659 3170 47.3 99.3 299,779 
11 399 .870 37.4659 3170 47 .3 99.3 299,779 
12 399 .870 37 .4659 3370 44.6 99.3 299,776 
13 399 .872 37 .4659 3270 45.9 99.3 299,779 
14 399 .870 37.4659 3270 45.9 99.3 299,777 
15 399 .870 37.4659 3270 45.9 99.3 299,778 
16 399 .870 37 .4659 3270 45.9 99.3 299,77 





*Cavity dismantled and contacting surfaces of end plate copper-plated. 


Since so much difficulty had been experienced with contact resistances, it 
was apparent that some rough estimate of their disturbing influence on Q, 
and on c, would be desirable. After the above determinations were com- 
pleted, the contacting surfaces of both end plates and cylindrical conducting 
surfaces were lightly amalgamated with mercury, and the cavity reassembled. 
The Q was then measured and was found to have increased to about 4400. 
Although the cavity dimensions had been altered slightly by the small fillet 
of mercury at the junction and spreading of the film of mercury over the interior 
of the cavity within about 2 mm. of the junction, six further measurements of 
c were made. These were all in the range 299,774-777 km./sec., in fair agree- 
ment with the previous set. However, although the results derived from the 
use of mercury contacts are inconclusive, the presence of residual contact 
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resistances in the first set of determinations is clearly demonstrated. Their 
effect on the measured value of ¢ is uncertain because the magnitude of the 
probable alteration of cavity dimensions when the mercury was applied to the 
junctions is indeterminate. Consequently, these measurements were con- 
sidered only qualitatively. 


(D) Errors 
Random 

Random errors entered the determinations in the measurement of cavity 
length and diameter constancy, the estimate of the center frequency of the 
resonance curve, the determination of the Q, the correction due to the refrac- 
tive index of air, and the temperature. Since the length measurement was 
carried out only to +23 yw, and since further it could be repeated to this accu- 
racy, it was judged to give rise to an error of about +2 km./sec. Errors due to 
inaccuracies of diameter measurement of the cylindrical conducting surfaces 
probably did not exceed +1 km./sec. 

Point scatter and small asymmetries of the resonant curve limited the 
accuracy of determination of the center frequency of any curve to about 
4 kc./sec. in the resonant frequency of about 400 Mc./sec. This led to an error 
of about +3 km./sec. in c. In Fig. 2 is shown a typical resonance curve pro- 
viding one of the 16 measurements. 

A random error in ¢ also entered due to errors in the determination of the 
cavity Q. The expression in the theory defining the resonant frequency of a 
coaxial cavity may be written 


c= 2 fo nL( (1+1/2Q). 


It is thus apparent that for small Q, the correction term involving 1/2Q 
becomes appreciable. In the present case, this term was in the vicinity of 
50 km./sec. Since Q was determined from the resonant curves by the expression 


Q = fo/(fe—fi) 


where f2 and f; are the half power frequencies, point scatter and non-linearity 
of the power detector led to an error in Q. From the experimental data this 
was estimated to be about +3%, so that the correction term to ¢ may have 
been in error by about +1.5 km./sec. 

A further error entered the determination of c because of inaccuracies in the 
determination of the refractive index of air, which depends on the air tempera- 
ture, pressure, and humidity. The first two variables could be measured to a 
sufficiently high degree of precision that errors from these were negligible. 
However, dependence of the refractive index of air on humidity is very marked, 
and uncertainty in this variable (measured by wet and dry bulb temperatures) 
was considered sufficient to lead to an error of about 1.5 km./sec. 

Errors in temperature measurement, which would affect cavity length, were 
estimated to result in an error in ¢c not more than +0.5 km./sec. 

The total maximum possible random error of a single determination thus 
may be as high as +9.5 km./sec. The square root of the sum of the squares 
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of these errors above is +4.3 km./sec. From Table I we can, assuming normal 
distribution, compute the probable error of a single observation as +2 km./sec., 
but this number is interesting only as an indication of consistency. 


2. Systematic 

The most serious systematic error, and the one whose elimination appeared 
the least likely, was that due to the contact resistances. It appeared from the 
measurements that the better the contact, the higher the value of c, but for 
values of Q above 3000 there did not seem to be any observable correspondence 
between Q and c, though perfect contacts might have increased the measured 
value of c still further. Tighter fitting of the end plates in order to attain better 
contacts was not permissible because the resulting distortion of end plates and 
cavity walls encountered when the plates were forced on would produce 
further errors. The tightness of fit used in the above determinations was such 
as to cause some distortion, especially on the inserts, but not sufficient to add 
seriously to the errors. 

Systematic errors also arose from cavity wall perturbations, probe holes, and 
probe penetration. The first was measured and found to be negligible although 
the diameter measurements were deemed to introduce the random error 
previously mentioned. The second was assumed to be small, since the ratio of 
probe hole area to cavity wall area was about 4 X 10-7. The last was effectively 
eliminated by retraction of the probes to a point beyond which further retrac- 
tion did not affect the resonant frequency. 

Another source of error, the magnitude of which was difficult to evaluate, 
was the change of surface conductivity due to work hardening caused by the 
cutting, grinding, and polishing operations. This would tend to increase the 
electrical length of the cavity, make our length determination too short, and 
thus result in a low value of c. This effect is mentioned by DuMond and Cohen 
(4) and if present could increase c by several km./sec. 


CONCLUSIONS 


It is apparent that the results of this research do not provide a definitive 
value of c. The principal source of uncertainty, which is inherent in the choice 
of a coaxial cavity excited in its TEM mode, lies in the fact that surface cur- 
rents cross boundaries between the ends and the walls of the cavity, and it 
would be extremely difficult to construct a coaxial cavity in such a manner 
that junctions would occur parallel to lines of current flow. The use of other 
types of cavities not subject to these conditions is more or less precluded by the 
large dimensions demanded at the low frequency employed. The wall conduc- 
tivity, not being perfectly uniform, will affect the correction term 1/2Q. This 
and other factors mentioned above lead the authors to believe that the value 
of 299,780 km./sec. is too low, containing an indeterminate systematic error. 
The estimated random errors remain superimposed. In our opinion these results 
point toward using the higher value of other investigators of 299,792. Experi- 
mental limitations of the accuracy attainable do not appear to warrant further 
work on coaxial cavities at this frequency. 
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AN APPARATUS FOR COMPARISON OF THERMOCOUPLES' 


By T. M. DAUPHINEE 


ABSTRACT 


This paper describes a semiautomatic apparatus for routine or precision com- 
parisons of thermocouples of the same type in the temperature range 0-1100°C. 
The couples being compared are welded together at the tips and placed in a tube 
furnace which is heated at rates varying from 10 to 100°C./min. Measurements 
of carefully annealed thermocouples show that in the temperature range 300- 
1100°C. platinum - platinum 10% rhodium thermocouple comparisons may be 
made to accuracies of +0.3 wv. (+0.03°C.) at heating velocities as great as 
15°C./min. while accuracies of +1.5 uv. at velocities of 100°C./min. are feasible. 
The furnace temperature is varied by means of a motor-driven variac with 
automatic reversal at peak temperature. In addition to this standard comparison 
procedure, provision is made for comparing corresponding elements of the 
couples, for suppression of all or part of the measured e.m.f., and for measuring 
the whole e.m.f. of a// couples when a comparison of different types is desired. 
The system can be adapted to XY recording with total e.m.f. plotted against 
e.m.f. differences. 


INTRODUCTION 

Calibration of thermocouples is ordinarily done in one of two ways: either 
by direct determination of the e.m.f. at thermometric fixed points, or by 
comparison of the couple with a standard one already calibrated. While the 
former method is precise, it is slow, and the thermocouples must be calibrated 
one at a time. The second method is potentially much quicker and quite accu- 
rate enough for most calibrations for industrial purposes. In addition it has 
wide possible applications in studies of the effect of heat treatment and ther- 
mocouple stability and reproducibility which require large numbers of high 
quality thermocouples of the same type to be treated. 

The usual procedure in comparing thermocouples (2) has been to weld them 
together at the tips (to ensure equality of temperature of the measuring 
junctions) and insert them in some sort of tube furnace which is then brought 
to an almost steady temperature state. The e.m.f.’s of the two couples are then 
read simultaneously, the process being repeated at a number of temperatures. 
The whole e.m.f. of each couple must be read to a large number of significant 
figures because the tips are electrically connected and the couples cannot be 
wired for difference measurements in the normal way. However, with thermo- 
couples of the same type it is possible to measure the e.m.f.’s between both sets 
of corresponding elements and from these the difference between the e.m.f.’s 
of the two thermocouples can be determined. This procedure has advantages 
since the differences between corresponding elements are small and change 
relatively slowly with temperature, so that a less precise knowledge of the total 
e.m.f. is required and higher heating velocities can be tolerated. Unfortunately, 
the measuring error will be doubled since two readings are required and the 
determinations must still be made simultaneously. 


1Manuscript received January 27, 19565. 
Contribution from the Division of Physics, National Research Council, Ottawa, Canada. 
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Another way of getting directly to the required difference at high heating 
velocities is by using the isolating potential comparator circuit (1) as shown in 
Fig. 1. A subtraction of the two thermocouple e.m.f.’s can be made by using a 


STANDARD , we caer ee TEST 


THERMOCOUPLE r--- ‘wee -5 THERMOCOUPLE 





wood 


Fic. 1. The basic thermocouple comparison circuit. 


double-pole, double-throw chopper (vibrator switch) and a condenser to 
transfer the e.m.f. of one thermocouple into series opposition with the e.m.f. of 
the other; the difference can be read with a potentiometer, using a galvanometer 
or other detector as indicator. If the potentiometer is of the self-balancing 
recorder type it will readjust itself until E, = E,+£,, where £, is the e.m.f. of 
the standard thermocouple, E,; is the e.m.f. of the test couple, and E, is the 
e.m.f. produced by the potentiometer and read on the recorder chart. The 
recorder reading is therefore a direct measure of the difference in e.m.f. of the 
two couples. The measurement of this difference is practically independent of 
rates of change of the whole e.m.f. provided the difference itself does not change 
too rapidly at the same time. As the diagram shows, two chopper and conden- 
ser units should be used at 180° phase difference to obtain the highest sensiti- 
vity. 

The whole e.m.f. of the standard couple may be measured to sufficient accu- 
racy for calibration purposes (approximately the same number of significant 
figures as is required of the difference) with relatively simple equipment. This 
method of comparison can be readily adapted to provide automatic calibration 
of a number of thermocouples simultaneously, since a maximum of three sig- 
nificant figures is sufficient for all required quantities and the number of func- 
tions to be recorded is small. 

The apparatus used in this laboratory was developed primarily for automatic 
comparisons of similar thermocouples in the 0-1100°C. temperature range. 
The system can compare and record the e.m.f. differences between one or more 
standard thermocouples and up to eight test thermocouples in a single tem- 
perature cycle from room temperature to 1100°C.-and back. Once started it 
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requires no further attention from the operator. Up to three couples can be 
fully analyzed, i.e. a comparison of the individual elements of the couples is 
also recorded. Provision has also been made for comparing thermocouples of 
different types (e.g. Pt- Pt, 10% Rh with Pt— Pt, 18% Rh or chromel- 
alumel) and for recording data in fixed point calibrations. 

THE THERMOCOUPLE CIRCUITS 
Comparison of Thermocouples of the Same Type 


The basic components of the thermocouple comparison circuit are illustrated 
in block form in Fig. 2. The hot junctions of a number of thermocouples are 
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Fic. 2. A block diagram of the actual thermocouple comparison circuit. 
welded together and heated in a nichrome tube furnace, the cold junctions 
being placed in an ice bath and attached to copper extension leads. The 
extension leads are taken to a terminal board which is in turn connected to a 
selector switch. This switch determines the quantity to be presented to the 
recorder, either directly or through the high sensitivity linear amplifier. 
Potentiometers (pot.) may be used for scale expanding purposes and an iso- 
lating potential comparator circuit associated with the switch allows differences 
to be presented when desired. 

The circuits used in presenting the various quantities to the recorder are 
shown in Fig. 3. When making difference measurements the functions are 
presented in sequences of three, the difference between one thermocouple and 
the standard couple (Fig. 3a), the whole e.m.f. of the standard (Fig. 35), and 
the difference between another thermocouple and the same standard. The 
switching is done approximately every ten seconds by the use of a solenoid 
drive on a Leeds and Northrup (31-3) 12-position selector switch modified 
in such a way that it will bridge between five different pairs of leads in any 
position.* The terminal board is laid out to facilitate cross connections where 
desired, and it is an easy matter to wire thermocouples so that they appear in 

*The modification consists of taking a five-pole selector of the usual type and replacing each of the 


five single tap (slip ring) sections by a 12-contact one. The switch brushes then join five different 
pairs of terminals on the switch at each switch position. Wiring ts shown in Fig. 4. 
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Fic. 3. Circuits used in the comparison of two thermocouples of the same type. In all cases 
S refers to the standard thermocouple and X to the thermocouple under test. The subscripts 
+ and — denote the positive and negative elements respectively. (a) Difference measurement. 
(b) Total e.m.f. measurement. (c) Comparison of corresponding elements (S and X) using the 
isolating potential circuit. (d) Comparison of corresponding elements by direct measurement. 


more than one position on the switch. For instance, a common practice is to 
show one group of three thermocouples four times per revolution. Poten- 
tiometers are not required for comparisons of identical thermocouples, but are 
needed for comparisons of, say, Pt - Pt, 10% Rh vs. Pt — Pt, 18% Rh. Two 
Rubicon type B potentiometers are permanently installed, another instrument 
being brought in when required, but for much of the work potentiometers are 
unnecessary and all the connections go directly to the amplifier. 

The amplifier is a Leeds and Northrup microvolt model giving 10 mv. output 
for 50, 100, 200, 500, 1000, and 2000 uv. input respectively on its various 
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Fic. 4. Wiring diagram showing the connections to the selector switch. 


ranges. It has 10 megohms effective d-c. input impedance and therefore does 
not load the thermocouple circuit appreciably. The record is taken by a Leeds 
and Northrup Speedomax single-point recorder with continuously variable 
sensitivity from 1 to 20 mv.* full scale, a small portable potentiometer being 
used to shift the zero to center scale for difference measurements. The chart 
normally runs at 1 ft./hr. with six seconds of pen record per point; a solenoid 
operated pen lifter causes a four second blank on the chart while switching takes 
place. 

The chopper unit is composed of four Leeds and Northrup 3338-6 choppers 
(single-pole double-throw, break before make) driven synchronously at about 
16 v., 50 c.p.s., by a common power supply and connected to give two double- 
pole units 180° out of phase with one another. Polystyrene or mylar condensers 
(about 100 uf.) are used; paper condensers cannot be used as they may intro- 
duce errors due to dielectric absorption. 


Comparison of Corresponding Elements 

The circuit represented by Fig. 3a can also be used to give the e.m.f. differ- 
ence between corresponding elements of the various couples. In this case the 
leads are connected on the terminal board as indicated by Fig. 3c. The whole 
e.m.f. (3b) and the difference between the two couples (3a) are also put on the 
chart, thus giving a complete record of the relationships between the couples. 
Three thermocouples can be compared in one run. 

Corresponding elements can also be compared by connecting them so that 
the whole e.m.f. of the standard goes to the recorder while the differences 
between corresponding elements are presented to the amplifier (3d). However, 
the difference between the thermocouple e.m.f.’s cannot be measured at the 
same time. 


Comparison of Thermocouples of Different Types 


The circuits described above are particularly useful for comparing thermo- 
couples that are of the same type or at least have similar sensitivities and 


*This range could be changed for greater convenience to, say, 2.5 to 50 mv. 
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e.m.f.’s. For thermocouples with widely differing sensitivities, e.g. Pt — Pt, 
10°0 Rh vs. chromel-alumel, total e.m.f. measurements are preferable and 
provision is made for easy conversion to this form of operation. The connec- 
tions are altered to permit total e.m.f. readings of all thermocouples through 
the three potentiometers and the amplifier. Two types of thermocouple may 
then be calibrated against a standard. For calibration purposes the poten- 
tiometers are changed by discrete steps of, say, 0.1 mv. and the recorder scale 
is used for interpolation. The same connections can be used for calibrations of 
couples at the various fixed points, up to three calibrations being made at one 


time. 


Other Forms of Operation 

Two Leeds and Northrup type HS galvanometers with lamps and scales 
have been installed and may be used with the potentiometers in the usual way. 
A completely separate set of leads from the terminal board is provided and 
neutral push-button type switches (0.01 uv. thermal e.m.f.) are used to select 
the thermocouples. 

The apparatus can also be used to record temperature differences with 
thermocouples, even in those cases where the sensing elements must be in 
direct electrical contact by being welded to the same piece of metal. The 
reading of the temperature difference is independent of temperature except 
as it affects the-thermocouple sensitivity.* 


Elimination of Errors 

Considerable care has been taken to prevent spurious e.m.f.’s from causing 
errors in the readings. The cold ends of all thermocouples are taken to the ice 
bath where they are joined to light leads of selected copper wire. These exten- 
sion leads, which are permanent connections in the case of our own couples, are 
joined to heavy duty wires, via copper terminal blocks, at the terminal board 
which is located near the furnace; lighter wire is used for the connections at the 
thermocouple switches. The switches and filter chokes (Fig. 5) are wholly 
enclosed in an insulated metal box to prevent temperature gradients from being 
built up in them. Low thermal solder is used for all soldered connections in the 
thermocouple and amplifier circuits and the wiring is arranged so that all 
soldered joints are in either the ice bath or the insulated box. The poten- 
tiometers and other circuit components are shielded from all sources of heat 
including the furnaces. When work on Pt — Pt, 10% Rh thermocouples requires 
maximum precision, the potentiometers are removed from the circuit entirely 
so that the small (up to 0.3 wv.) thermal e.m.f.’s normally found in them will 
not cause errors. As a precaution against inhomogeneities arising from bending 
the thermocouples, all couples are mounted on a small wheeled carrier, which 
also carries the ice bath for the cold junctions. This carrier is lined up with the 

*Note added in proof.—Since submission of this paper it has been found that because of the 
speed with which thermocouples can be handled it 1s never necessary to compare more than six at 
one time (two standards and four unknowns). Consequently the terminal board has been reduced 
to the size shown in Fig. 7. This has enabled us to use one of the push-button switches (which are 


located below the galvanometer scales) to select the type of operation desired. Cross connections on 
the terminal board have been entirely eliminated. 
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Fic. 5. The filter system. 


furnace so that the thermocouples can be inserted or withdrawn without any 
flexing whatsoever. 

All parts of the measuring equipment are electrically shielded and a common 
ground system connects them with the furnaces. Heavy filtering (Fig. 5) is 
used to isolate the choppers from a-c. pickup arising in the furnaces and the 
power transformers. The furnaces are powered by a-c. through line isolating 
transformers and have been provided with optional hum balancing ground 
connections. These do not appear to be required but are useful for checking 
purposes. 

POWER AND CONTROL CIRCUITS 


The Furnace and Furnace Power 

Comparisons are made in a standard tube furnace of the type shown in 
Fig. 6. The 15 kva. motor-driven variable autotransformer is supplied from the 
220 v. a-c. line. Two gear ratios, with a factor of difference of about 100, are 
provided and the speed of the motor itself can be varied. Push buttons control 
the raising and lowering of voltage, which can be either continuous or for the 
duration of the push. In addition there is provision for automatically stopping 
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Fic. 6. The tube furnace and power circuit. 
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the motor at full and zero voltage, for reversing the motor when the recorder 
reaches a preset value, and for giving a warning signal at the maximum of the 
recorder scale. The motor automatically turns the transformer to zero at full 
speed when the apparatus is shut off, and immediately after a power inter- 
ruption. Interlocks are provided so that the furnace main switch opens auto- 
matically and stays open if a-c., d-c., or water supply should fail, and at 
night a timer switch shuts off both a-c. and d-c. power. As a further safety 
precaution against burnout, the timer and main a-c. switch are of the manual 
reset type. The tubes carrying the cooling water from one end of the furnace 
tube to the other pass through the core of the driving transformer alongside 
the single turn secondary and therefore do not require an insulated center 
section. A ground balancing potentiometer (25 ohms) is mounted on the side 
of the furnace. It is connected to the ends of the nichrome tube and has its 
center tap grounded. 


Switch Drive and Pen Lift Mechanism 

The selector switch is driven in 1/12 revolution steps by a ratchet and pawl 
mechanism which is gravity operated with a solenoid return. Switching time 
may be adjusted from 3 to 20 sec. by varying the constants of the RC timing 
network employed, and an adjustable delay is provided between lift and fall 
of the mechanism so that the pen lifter will have time to operate. 





Fic. 7. The control panel and measuring equipment. 
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Fic. 8. A complete comparison of three thermocouples from different manufacturers taken at slow speed (1 
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Layout and Controls 

The apparatus is laid out in three areas. The power units (i.e. the variable 
voltage transformers with the motor drive, the interlock safety switches, and 
the main switch) are in one area, the furnaces and terminal board are in another, 
and the measuring circuit and control panel are in the third. A 6 to 8 ft. spacing 
between the different groups of components minimizes the effects of stray 
fields from the transformers and furnaces. Fig. 7 shows a photograph of the 
tube furnace and terminal board and of the measuring and control section. 


RESULTS 


The apparatus has so far been used chiefly for intercomparisons of standard 
Pt — Pt, 10% Rh thermocouples in the temperature range where these couples 
are used to define the International Temperature Scale. Heating rates have 
varied from 5 up to 200°C./min. and a large number of complete heating and 
cooling curves have been obtained. A complete comparison (the differences 
between corresponding elements also shown) of three of these thermocouples, 
all from different manufacturers, is shown in Fig. 8. A comparison of results 
for different heating rates is shown in Fig. 9. For carefully annealed and 
treated thermocouples the error arising from the fact that the temperature is 
changing is normally less than 0.2 wv. at 15°C./min. and is rarely as much as 
1 wv. at 100°C./min., except at very low temperatures where the junctions lag 
a considerable amount behind the furnace. Averaging the up and down runs 
appears to result in an over-all accuracy for the measuring equipment of about 
+0.3 wv. and 1.5 pv. respectively. It seems probable that a large part of the 
residual errors is contributed by imperfections in the thermocouples themselves 
or by gradients in the furnace. 

The comparisons obtained with this apparatus have an accuracy comparable 
with the best that can be attained by the two potentiometer method and are 
obtained with considerably less expenditure of time and effort. They have the 
further advantage that a permanent record of the calibration data is obtained 
which can be very useful in checking for reading errors. The accuracies 
required of calibrations for most industrial purposes are such that very con- 
siderable heating rates can be used, and the number of couples that can be 
processed is limited largely by the time required in preparing them for the run 
and in reading the charts. The use of a function plotter (X Y recorder) could 
introduce a further saving of time in chart reading and plotting. 
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AN INVESTIGATION OF THE NUCLEAR RESONANCE 
ABSORPTION SPECTRUM OF Al” IN A SINGLE 
CRYSTAL OF EUCLASE! 


By R. G. EApDEs? 
ABSTRACT 


The nuclear magnetic resonance spectrum of Al? in a single crystal of 
HBeAISiO; (euclase) in an external magnetic field of 7000 gauss is described. 
The Al?’ nuclei are found to occupy two kinds of lattice sites which differ from 
each other only in the orientation of the principal axes of the electric field 
gradient tensors whose eigenvalues are all distinct, but identical at both sites. The 
absolute value of the quadrupole coupling constant for both groups of Al” nuclei 
is found to be 5173+10 kc./sec. and the axial asymmetry parameter of the field 
gradient tensors at both the lattice positions is 0.698+0.009. The two sets of 
principal axes are reflections of each other in the (010) plane which is consistent 
with the symmetry of the proposed crystal structure. The theoretical dependence 
of the spectrum line frequencies on the applied magnetic field is given for one 
group of AJ?’ nuclei and for one orientation of the crystal in fields up to 2000 
gauss; the feasibility of checking this experimentally is discussed. 

1. INTRODUCTION 

Volkoff and Lamarche (3, 7) have investigated theoretically the nuclear 
magnetic resonance spectrum of the Al?’ nucleus (J = 5/2) placed in an 
electric field gradient which is not axially symmetric. Their results were stated 
explicitly for the case of spodumene where the large asymmetry parameter 
n = 0.95 (0 < n < 1) of the electric field gradient enhances the departure of 
the line frequencies from a linear dependence on the applied magnetic field 
H for fields less than about 800 gauss. The relatively low quadruple coupling 
constant of Al?? in spodumene, eQ¢,,/h = 2.960 Mc./sec., results in a pure 
quadrupole spectrum (H» = 0) of two lines at 789 and 758 kc./sec. and a weaker 
line at the sum of these two frequencies. This places the whole of the spectro- 
scopically interesting region for which perturbation calculations are no longer 
valid in a frequency range below 1 Mc./sec. So far, experiments made to 
observe spectra in this frequency region have not been successful and the 
predictions of (3) and (7) have yet to be checked. Accordingly a search has 
been made for a crystal in which the Al?’ nuclei are at lattice sites where the 
field gradient tensor lacks axial symmetry and in which the quadrupole coup- 
ling constant is somewhat greater than in spodumene. 

A number of aluminum silicates have been examined and the results for 
one of these, euclase, are given in detail here. The theoretical treatment given 
by Volkoff in (6) is used to analyze the experimental data and the possibility 
of using this crystal to check the Volkoff-Lamarche predictions is discussed. 


2. THE EUCLASE CRYSTAL 
Euclase (HBeAISiO;) is a monoclinic mineral crystal with space group 
C3,. Biscoe and Warren (2) give a projection of the structure on the (a, c) 
plane. The parameters of the unit cell are given by them as: 
a=462A, b= 14.24A, c=4.75A, B = 79° 44’. 


'Manuscript received March 29, 1955. 

Contribution from the Department of Physics, University of British Columbia, Vancouver, 
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The structure they proposed is characterized by independent SiO, groups and 
is considered to belong to the orthosilicate group. Each aluminum atom is 
surrounded by five oxygen atoms and one OH group which together form dis- 
torted octahedral groups. These groups form zigzag chains which run parallel 
to the c axis of the crystal. This proposed structure is supported by the exis- 
tence of the perfect cleavage parallel to the plane (010). 

The aluminum nuclei do not lie on symmetry axes of the crystal as in the 
case of spodumene; neither of the two experimental crystal rotations could 
therefore be made about a known principal axis of the field gradient tensor. 
The two rotational axes were instead chosen with respect to the cleavage 
parallel to the plane (010). 

Two single crystal specimens of euclase were used: both were pieces of larger 
crystals and on each the (010) cleavage was excellent. A slightly tinted crystal 
about 12 X 8 X 6 mm. in size was used for the Y or 0 axis rotation, the crystal 
being mounted so that it could be rotated about an axis perpendicular to the 
cleavage face. This crystal was also used for the initial observations on the X 
axis rotation. Later, a larger clear specimen was obtained and this was used 
for the final observations on the X rotation. Examination of this specimen 
showed that the X axis was, in fact, parallel to the c axis of the crystal. The 
sets of data obtained from these two crystals agreed within the experimental 
uncertainty but the results quoted for the X rotation are those obtained with 
the clear crystal since this both increased the available signal to noise ratio 
and enabled the X axis to be more accurately defined. By using carefully 
machined surfaces on the crystal mounts and suitable jigs the two rotation 
axes were made mutually perpendicular to within 0.1 degrees. 


3. APPARATUS 


The recording nuclear magnetic resonance spectrometer constructed was 
similar in basic design to many now in use. (See for example Pound and 
Knight (5).) The external magnetic field was supplied by a large permanent 
magnet having 8 in. diameter pole caps. These caps were in the form of trun- 
cated cones and follow the design data published by Andrew and Rushworth 
(1). The field gap was 2 in. and in the best region the homogeneity was better 
than 0.2 gauss in 7000 gauss over a 2 cc. volume. 

The oscillating detector design followed closely that developed by Watkins 
(9). To achieve optimum in signal to noise and to permit an accurate angular 
measuring device to be incorporated, the detector was built on a heavy brass 
base-plate. This carried the coaxial arm leading down to the coil assembly in 
the field. The radio-frequency coil was wound directly on the crystal which was 
fixed with respect to the arm. This procedure ensured the best possible filling 
factor. A medium power eyepiece and mirror system was used to read the 4 in. 
diameter 360° scale mounted on the coaxial arm, which, for rigidity, was made 
of stainless steel tubing. 

The base-plate was provided with levelling screws which allowed the axis of 
rotation to be made perpendicular to the magnetic field. This was first done 
approximately using spirit levels and then more accurately by observing the 
spectrum at 180° points: frequency measurements on the component lines 
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should agree within the other limitations of accuracy. This adjustment was 
made using the Y rotation for which the spectrum has only half the number 
of lines observed for the X rotation and hence approximately twice the inten- 
sity. The small residual error from this source was then further reduced by 
averaging observations taken at 180° points. The accuracy of individual 
angular settings was +0.1 degrees. 

Frequency measurements were made with a Lampkin Frequency Meter, 
Type 105B, and an Eddystone communication receiver. This combination 
allowed the frequencies of individual lines to be measured from the recorder 
chart with an accuracy of +0.15 kc./sec. for the strongest lines and +0.50 
kc./sec. for the weakest lines. Variations of the magnetic field of the permanent 
magnet with temperature were significant and the Al?’ signal from a solution 
of AICI; was used to calibrate the magnet. The temperature coefficient was 
equivalent to a decrease in the Larmor frequency of the Al?’ nucleus of 1.28 
kc./sec./°C. rise in temperature. 


4. RESULTS 


In (6) Volkoff extends the original theoretical investigations by Pound (4) 
and Volkoff et al. (8) and gives explicit formulae for the dependence of the 
absorption frequencies on the angle of rotation of the crystal about any 
arbitrary axis perpendicular to the applied magnetic field. These formulae, 
applicable to crystals in which non-axially symmetric crystalline electric field 
gradients exist at the sites of the nuclei being investigated, are valid to second 
order terms in the quadrupole coupling constant eQ¢,,/h. In particular he 
shows that when the innermost satellite lines (m = 3/2, for I = 5/2) are 
visible and second order nuclear quadrupole effects are present, rotations 
about two arbitrary but mutually perpendicular axes yield complete informa- 
tion concerning the field gradient tensor at the site of the nucleus and also 
give the absolute value of the quadrupole coupling constant. Two such rota- 
tions were made. 


(a) Data from the X Rotation 


When the crystal was mounted with the long axis of the cleavage plane 
parallel to the axis of rotation the spectrum consisted of 10 lines of markedly 
different intensities with the weakest lines more than 2 Mc./sec. apart at 
particular values of 0x. The intensities of the lines were functions of the angle 
of rotation 6x. Since the maximum number of lines expected is 2/, or five for 
I = 5/2, it is evident that the Al?’ nuclei in euclase occupy two non-equivalent 
lattice sites. 

The frequencies of the four innermost satellite lines, called Ax, By, Cx, and 
Dx, are shown as functions of 6x in Fig. 1. Each plotted point represents the 
mean of at least four frequency determinations made from selected recordings. 
The points common to Ax and Cx, and By and Dx, respectively, were frequen- 
cies recorded at exact overlap; this point could be located to within +0.1 
degrees by observing the change in intensity of the combined lines. The fre- 
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quency differences Ey = (v4—vp)x and Fy = (ve—vg)x are plotted in Fig. 2 
and have the form 
axyt+by cos 26x + Cx sin 20x 


corresponding to equation [21] of Ref. (6). 
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Fic. 1. X rotation. Frequencies of the four innermost satellite lines of the Al?’ spectrum 
plotted as functions of the angle of rotation 6y. 
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Fic. 2. X rotation. Frequency differences (v4 —vp)x and (vc —vg)x. The open circles are the 
experimental points and the smooth curves are the functions 
Ex = 298.9—198.7 cos 26x +734.8 sin 26x, 


Fx = 298.9 — 198.0 cos 26x —731.5 sin 26x 
fitted to these points. 
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Fic. 3. X rotation. Perturbation shift of the two central components from the Larmor 
frequency. The open circles are the experimental points and the smooth curves are the functions 


Gy = —27.27+49.97 cos 20x —28.40 sin 20x —15.73 cos 40x —8.83 sin 46x, 
Hy = —27.15+50.18 cos 26x +28.72 sin 26x — 15.85 cos 40x +8.52 sin 40x 
fitted to these points. 


The data for the split central line (m = 3) are given in Fig. 3 where the 
perturbation shift from the Larmor frequency is plotted as a function of 4x. 
Second order effects are present and Gy and Hx have the general form of 
equation [23] of Ref. (6), i.e. 

Nx+px cos 26x+rx sin 20x+Ux cos 40x+vx sin 46x. 
Fourier analyses of Ex, Fy, Gy, and Hx gave the constants listed in Table I 
where all the symbols are as defined in Ref. (6). The error limits quoted in the 


TABLE I 


X ROTATION DATA: VALUES OF CONSTANTS IN KC./SEC. REQUIRED TO FIT 
CURVES TO EXPERIMENTAL POINTS 








Group 1 Group 2 
nuclei nuclei 
Center line Gy Hx 
Satellite lines Ay and Dx By and Cy 
ax 298 .9 298.9 +4.0* 
bx —198.7 —198.0 +5.3 
Cx 734.8 —731.5 44.7 
nx —27 .27 —27 .15+0.25 
bx 49.97 50.18+0.33 
rx — 28.40 28 .72+0.30 
ux —15.73 —15.85+0 .33 
vx —8.83 8.52+0.30 


*Calculated uncertainties are the same for both groups. 


data for both the X and Y rotations are based on a total initial uncertainty 
of 0.5 ke./sec. and 2.5 kc./sec. in the determination of the frequency of a 
central line and satellite line respectively. These estimates of uncertainty 
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include contributions from (a) the measurement of the angle of rotation, 
(6) the determination of the frequency, and (c) the value of the magnetic 
field (temperature effect). In the numerical calculations the uncertainties have 
been treated as probable errors. 

Using the Fourier amplitudes mx, px, rx, Ux, and vy obtained from the 
analyses of Gy and Hx, values of by and cy were calculated (Ref. (6), equation 
[24]) and compared with those obtained directly from the satellite pairs Ex 
and Fy. In this way the central line Gy was associated with satellites Ay and 
Dy, and Hy with By and Cx. 


TABLE II 


X ROTATION DATA: EIGENVALUES IN KC./SEC. AND ASYMMETRY PARAMETERS 
OF THE QUADRUPOLE COUPLING TENSOR Kq,; AT THE Al?’ SITES 





Group 1 Group 2 
Eigenvalue nuclei nuclei 
Kédz:z —234.4+8.0 —224.8+7.4 
Rides —1318.7+44.3 —1324.744.5 
K¢:: 1553 .0+3.5 1550 .0+3.0 
n 0.699 +0 .007 0.704+0 .008 





The first three equations of [24], Ref. (6), reduce to 
Cy*+cz* = 4.5ay?—1.75R? —Qvonx, 
Cy’ +z? = axbyt+3vopx, 
CyCz = 1.dv0rxy+0.5axcx. 

Using these relations, the Fourier amplitudes cy and cz were then calculated 
in terms of ax, by, Cx, Nx, Px, and rx. For the two groups of nuclei these were: 
cy: = +(1102.3+5.0), Cy2 = +(1103.6+5.0) kc./sec.; 

C21 = +( 203.7+3.7), Cz2 = +(193.943.9) kc./sec. 


The ambiguity in sign is removed by the Y rotation data which show that cy: 
and cy: are equal and both negative. The quadrupole coupling tensor K¢;; was 
then formed and diagonalized as discussed by Volkoff et al. in Ref. (8); the 
principal axes xyz were chosen so that ¢,, has the smallest and ¢,, the largest 
absolute value. » = (¢;:—¢,,)/¢2: is the asymmetry parameter and is thus 
restricted to the range zero to unity. The eigenvalues of the tensors for the 
two groups of nuclei are given in Table II along with the asymmetry parameter 
n. The process of diagonalization also gives, apart from an ambiguity in com- 
mon sign, the direction cosines of the principal axes xyz with respect to the 
set X YZ. The corresponding angles are given in Table IV along with those 
obtained from the Y rotation data. 


(b) Data from Y Rotation 


When the crystal was mounted with the axis of rotation perpendicular to 
the cleavage face the spectrum consisted of five lines. Again the intensities 








292 CANADIAN JOURNAL OF PHYSICS. VOL. 33 





Oo 20 40 60 60 100 120 10 160 180 200 
6y 
DEGREES 


Fic. 4. Y rotation. Frequencies of the two innermost satellite lines of the Al?’ spectrum 
plotted as functions of the angle of rotation 6y. 
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Fic. 5. Y rotation. Frequency difference (vg—va)y. The open circles are the experimental 
points and the smooth curve is the function 


Cy = —59.9+551.0 cos 26y — 1086.1 sin 26y 
fitted to these points. 


were functions of the angle of rotation 6y. In addition, the satellite lines were 
broader at positions where the shift from the unperturbed Larmor frequency 
was a maximum. The proposed crystal structure indicates that exact super- 
position of the pairs of lines from the two groups of nuclei only occurs with the 
b axis perpendicular to the magnetic field Hy and the broadening is therefore 
likely to be due to a small error in alignment. 

For the X rotation, 6x is zero when the Y axis is parallel to Ho, and the other 
crossover of lines Ax and Cx, or By and Dx, in Fig. 1 defines the position 
where the XY plane is perpendicular to Hy; for the Y rotation, this latter 
position defines 6y = 0. The crystal was mounted with the Y axis accurately 
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Fic. 6. Y rotation. Perturbation shift of the central component from the Larmor frequency. 
The open circles are the experimental points and the smooth curve is the function 
Dy = —44.63 — 19.45 cos 26y+11.04 sin 26y — 27.54 cos 40y —36.82 sin 46y 


fitted to these points. 


perpendicular to Hp and with 6y approximately zero. The complete data for the 
rotation were then taken and the relation 


(vay —¥dx) 6y=90° = (Vey —VBy) 6y=90° = (Vay —Vay) 6y=0° 


used to correct the zero of the 6y scale. 

The frequencies of the two innermost satellite lines A» and By are shown as 
functions of 6y in Fig. 4. The difference (vgy—va,) is plotted in Fig. 5 where 
the smooth curve drawn has been synthesized from the Fourier components 
obtained by analyzing the experimental data. The perturbation shift of the 
central line is shown in Fig. 6 where a marked 46, component is evident. 

Analysis of the experimental data follows that outlined for the X rotation 
and is summarized in Table III which gives the Fourier amplitudes of the 
curves fitted to the experimental points together with the calculated constants 


TABLE III 


Y ROTATION DATA: VALUES OF CONSTANTS IN KC./SEC. REQUIRED TO FIT CURVES 
TO EXPERIMENTAL POINTS, AND THE EIGENVALUES AND ASYMMETRY PARAMETER 
OF THE QUADRUPOLE COUPLING TENSOR 











Fourier constants Derived constants 
Satellite data || 
ay —59.9+4.0 | Cz +219.0+ 8.3 
by 551.0+5.3 cx +742.1412.8 
cy —1086.1+5.2 {| 
| 
Central line data 
Korz —240.2+ 9.3 
ny —44.63+40.25 Kéyy —1311.0+ 5.0 
py —19.45+0.33 Kez: 1552.0+ 4.0 
ry 11 .04+0.30 
uy —27 .54+0.33 
vy —36.82+0.30 ” 0.692+ 0.012 











294 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


TABLE IV 


ANGLES BETWEEN THE PRINCIPAL AXES xyz AND X YZ; THE MEAN 
PROBABLE ERROR IN THE ANGLES IS +45’ 


Angle 





between From Group 1 Group 2 
axes rotation nuclei nuclei 
xX X 112° 45’ 113° 5’ 

ye 112° 46’ 112° 46’ 
6, X 150° 56’ 29° 0’ 
Y 150° 20’ 29° 40’ 
Ment X 106° 52’ 106° 35’ 
y 107° 55’ 107° 55’ 
y, X X 143° 54’ 143° 52’ 
ig 144° 42’ 144° 42’ 
y) F X 79° 57’ 100° 17’ 
Y 81° 59’ gs° 1’ 
y, 2 X 55° 51’ 55° 48’ 
¥ 55° 55’ 55° 55’ 
Bk X 63° 54’ 63° 52’ 
y 63° 43’ 63° 43’ 
a X 116° 58’ 63° 12’ 
y 118° 23’ 61°37’ 
Coe xX 39° 11’ a0°° bY 
y 40° 31’ 40° 31’ 


cz and cx, the eigenvalues of the quadrupole coupling tensor A¢;,, and the 
asymmetry parameter. The ambiguity in the signs of cz and cx is significant 
here and gives the two sets of principal axes which correspond to those obtained 
from the X rotation data. Table IV gives for both rotations the angles be- 
tween the principal axes xyz and the set XYZ. 


5. DISCUSSION 

The only X-ray data available for euclase are those quoted in (2). Biscoe 
and Warren proposed a structure which was consistent with their experimental 
data and in agreement with empirical rules for silicate structures. But they did 
not claim a unique solution. The X (c axis) and Y (6 axis) rotation spectra 
obtained in this investigation confirm the symmetry elements of the proposed 
structure. With the a or c axis perpendicular to the applied magnetic field 10 
lines are expected in the Al?’ spectra as found for the X rotation. For the 6 or 
symmetry axis perpendicular to the field only five lines are expected and this 
too is confirmed by the Y rotation data. 

No correlation is apparent between the directions of the principal axes of 
the electric field gradient tensor (see Table IV) at the sites of the aluminum 
nuclei and the positions of any of the nearest neighbors. The environment of 
the aluminum nuclei has no great symmetry and at best the bond directions 
for the five oxygen atoms and the one OH group can only be found approxi- 
mately, since the a and c coordinates are stated by Biscoe and Warren to be 
only approximate. Further there is no way of allocating the correct tensor 


of the pair to the appropriate group of aluminum nuclei. 
The quadrupole coupling constant of Al? in euclase is 5.173+0.010 Mc./sec. 
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and is therefore appreciably higher than the value for spodumene although 
of the same order of magnitude. 

Lamarche and Volkoff have shown in (3) that the secular determinant of the 
Hamiltonian representing a magnetic nucleus (J = 5/2) in an applied mag- 


5.0 


3.0 


MC/sec 


°o ' 2 R 3 4 
° 466 932 Ho 1398 1864 GAUSS 


_ Fic. 7. Theoretical spectrum line frequencies in Mc./sec. for low values of the magnetic 
field Ho (or of R = 4pHo/eQ¢::). 


netic field Hy» and subjected to a crystalline electric field gradient can be re- 
duced to two 3 X 3 determinants for the special case in which Hp coincides 
with one of the directions of the principal axes of the electric field gradient 
tensor. The solutions of these determinants give directly the energy levels. 
The transition frequencies and their probabilities are then readily calculated. 
Using the notation of (3) the nine spectrum line frequencies of the Al?’ reso- 
nance have been calculated for the special case of the z principal axis along Hp. 
These values are therefore only applicable to the one group of Al?’ nuclei 
which satisfy this condition. For the other group none of the principal axes 
will be along the field direction and these nuclei will give rise to a maximum of 
15 lines of reduced intensity since transitions between any of the six energy 
levels are possible. In Fig. 7 the transition frequencies are plotted as functions 
of the applied magnetic field for values up to 2000 gauss. Table V gives actual line 
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frequencies for specific values of the dimensionless parameter R = 4uH)/eQ¢.:. 
The pure quadrupole resonance (H» = 0) consists of two strong lines at 1.11 
and 1.44 Mc./sec. and a weak line at 2.55 Mc./sec. 

An attempt is being made in this laboratory to verify these predictions. The 
accuracy of the entries in Table V is certainly better than +10 kc./sec. and a 


TABLE V 


TRANSITION FREQUENCIES IN Mc./SEC. FOR ONE GROUP OF AI?” NUCLEI IN 
EUCLASE AS A FUNCTION OF R: LINES » TO v5 ARE THE ZEEMAN LEVEL 
TRANSITIONS FOR HIGH VALUES OF THE APPLIED MAGNETIC FIELD Ho 
(R = 4uHo/eQd::) 





R vy ve v3 V4 V5 v6 v7 v3 V9 
0 0 1.44 0 beet 0 1.44 cSt 2.55 2.55 
0.25 0.64 1.26 0.31 1.25 0.06 1.59 1.00 2.21 2.91 
0.5 2d 1.08 0.61 1.42 0.11 E73 0.93 1.89 3.27 
0.75 1.83 0.938 0.80 1.56 0.13 1.97 0.89 1.69 3.65 
1.0 2.01 LI 0:57 Lae 0.11 2:56 0.91 1.92 4.04 
1.5 2.29 1.46 0.07 1.03 0.08 3.82 1.08 2.56 4.82 
2.0 230 eet 0.70 0.79 0.15 5.03 1.34 3.26 5.67 
3.0 3.08 2.32 1.54 0.90 0.06 6.95 2.49 4.75 7.90 
4.0 3.59 2.85 2.09 1.35 0:53 8.55 3.96 6.30 10.39 


search need only be made over this range. Pure quadrupole resonances have 
already been observed in other substances around | Mc./sec. and a check can 
therefore be made immediately at Hy) = 0. 

Previous experimental work with spodumene indicates that about | Mc./sec. 
is the lower limit of operation of the spectrometer in use, but further develop- 
ment may well reduce this considerably. Fig. 7 shows, however, that there are 
interesting features in the 1 to 3 Mc./sec. range for this crystal. There is 
therefore every expectation of checking the field-frequency dependence of 
some of the nine lines. 

Volkoff and Lamarche have extended the calculations of (3) and give in (7) 
the relative intensities of the spectrum lines for Al*? in spodumene for the 
radio-frequency magnetic field H, linearly polarized along the x principal axis of 
the electric field gradient tensor and the static magnetic field Hp along the z 
principal axis of this tensor. The interesting feature ot these results (which will 
be essentially the same for Al*’ in euclase) is that for small values of Ho both 
circularly polarized components of H, are effective in causing transitions. This 
results in an enhancement of the intensity of certain lines at the expense of 
others (Fig. 1, Ref. 7). It will clearly be best to look for the strongest lines 
first, two of the most interesting features being the marked increase in intensity 
of the line designated v2 around 1 Mc./sec. and the appearance of the line v7 
(not observed at high fields) with reasonable intensity around 2 Mc./sec. 
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MASSES OF NUCLEI IN ISOBARIC SEQUENCES! 


By KaAILASH KUMAR AND M. A. PRESTON 


ABSTRACT 


The semiempirical mass formula is known to give a rough fit to the general 
trend of atomic masses. A detailed correlation to within 0.5 Mev. or less is 
attempted between observed atomic mass-differences and the predictions of a 
semiempirical liquid drop formula generalized to include terms which depend on 
angular momentum in the way suggested by the nuclear shell model. For fixed A, 
it is found that, allowing for these ‘‘spin terms’’, the Z-dependence is parabolic. 
The variations of the different parameters with A and with shell changes are 
studied. It is found that, when close agreement with experimental values of mass 
is required, certain ‘‘constants”’ of the liquid drop formula are in fact quite de- 
pendent on A. This is interpreted to mean that the A- or Z-dependence of the 
formula is incomplete. Nevertheless, an empirical method of extrapolating un- 
known masses is suggested. 


INTRODUCTION 


A mass formula should be able to give (1) the 6-decay energies, (2) the 
atomic-numbers of 8-stable elements, and (3) the masses of the atoms. Various 
studies have been made with the liquid drop formula, and some sets of pa- 
rameters have been obtained* which give good agreement with the general trend 
of the masses. However, marked deviations from any of these formulae, of 
the order of several milli-mass-units, are known to occur at shell regions. 
Similarly, the line of B-stability obtained from such a formula is known to be 
displaced about a charge unit away from the observed values at certain places. 
On the other hand, there exist studies in B-energetics (4) which give a fairly 
good systematization of energies and stability, but the correlation of these 
with mass has not been investigated. Shell effects have also received attention 
and some ingenious formulae (7) have been put forward. These are, however, 
hard to interpret and even then are not enough to explain all variations. 
Presumably because of lack of data on 8-decays, no attempt seems to have 
been made to study all three aspects simultaneously in detail. 

In this work we have tried to correlate the information from these three 
sources. In Part I we use 6-decay data to evaluate the spin and pairing terms 
and the “parabolic’’ constants (B4 and Z,4). Allowance is made for large spin 
terms that arise in transitions across shells. The “‘parabolic’’ constants then 
obtained are independent of last particle or last pair effects and should be 
smooth functions of A. In Part II we try to fit the constants obtained in 
Part I and the absolute experimental masses within the framework of the 
liquid drop formula. In order to achieve exact agreement with all three criteria 
we find that at least two ‘‘constants”’ of the liquid drop formula must vary with 
mass-number. In other words, the dependence of the liquid drop formula on 
masses is inadequate to satisfy experimental requirements. Finally, by giving 

1Manuscript received February 10, 1955. 

Contribution from Physics Department, Hamilton College, McMaster University, Hamilton, 
Ontario. The work was partially supported by a grant from the National Research Council of 


Canada. 
* Most recent of these is Reference 6. Reference 8 gives the results of previous works. 
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up some of the liquid drop relationships between the constants, an empirical 
method of computing masses is suggested. 
I. RELATIVE MASSES OF ISOBARS 


1. Notation 
In the liquid drop or statistical model of the nucleus the variation in the 

masses of isobars is attributed to differences in coulomb energy, symmetry 
energy, and odd-even pairing energy effects. The first two of these are quad- 
ratic in Z, the charge number. Hence, with a slight generalization to include 
shell structure effects, we can write the mass of a nucleus A,Z as 
[1.1] M(A,Z) = M(A,Zo) + 2Bal(Z — 24)? — (Zo — Za)?] 

— D3 (2; — Zy)bs — Ly (Ny — Nyo)n, 

— 2005 (Ay — Ago)wy — 2s (Ay — Ajo); 


Here Z, is the charge of the lightest @-stable nucleus of mass-number A and 
the other symbols modified by subscript 9 refer to the values of the corre- 
sponding quantities in this stable nucleus. B, and Z,4 define the quadratic 
dependence on Z, B, determining the width of the Z-parabola and Z, the 
position of its minimum. The forms of B, and Z, in the semiempirical formula 
are discussed later, but we do not insist on this theory in this section. The 
summation terms in [1.1] take note of the less regular properties of specific 
nuclei not allowed for in the coulomb and symmetry energy terms. It is 
assumed that the binding energy and pairing energy can be considered con- 
stants within each nucleon subshell. Then the summation index 7 refers to a 
subshell and the symbols in [1.1] have the following definitions. Z,(N) is the 
number of protons (neutrons) in the j* proton (neutron) subshell of nucleus 
(A,Z); p; and m, are the corresponding binding energies per particle in the 
self-consistent potential of the shell theory, excluding the coulomb, symmetry, 
surface tension, and volume effects included in the semiempirical formula. 
In other words, p, represents the difference between the binding energy for a 
proton in the 7 subshell and the over-all average binding energy per particle. 
However, since our consideration involves only differences of pj's and n,j’'s, we 
shall refer to these loosely as the proton and neutron binding energies. \,; (A,) 
is the number of proton (neutron) pairs in the 7" subshell and 7, (v,) is the 
corresponding proton (neutron) pairing energy per nucleon. 

In accordance with normal usage, we shall refer to these summation terms 
as “spin” terms. 

In the following we assume simple transitions in which only the outermost 
nucleons are involved. In negative 8--decay of an odd-A, even-Z nucleus an 
even number of protons increases to an odd number, while an odd number of 
neutrons decreases to an even number. Thus, there is no change in the number 
of pairs and hence no contribution from the last two terms of equation [1.1]. 


Thus, the energy of transition is 
[1.2a] O% = LOZ+Ds—m, A odd, Z even, 


[1.3] where LQ¢ = B,(Z,4—Z—0.5). 
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In the Bt-decay* of odd-A, odd-Z nuclei, there is again no change in the 
number of pairs. For 8+-decay, we let Z denote the charge of the daughter 
nucleus and define Q% by QZ = —Q*(Z+1-—2Z). Then equation [1.2a] 
holds for these transitions also. 

In the B--decay of an odd-A, odd-Z nucleus or 8+-decay of an odd-A, even-Z 
nucleus, each nucleus involved has either a neutron or a proton pair more 
than the other. These transitions are represented by the equation 


[1.25] Of = LOZ+),—m, +27 j;—2%, A odd, Z odd. 
Similarly, 

[1.2c] Ot = LOZ+b;—m—2%, A even, Z even, 
[1.2d] 02 = LQZ+pj—m +27, A even, Z odd. 


These equations are valid only for transitions in which shells are filled regu- 
larly. For a general case, one has to work with equation [1.1], actually counting 
the number of pairs and particles. 

2. Computations and Treatment of Data 
For a fixed A, let us plot QZ against Z. This means that we plot B--decay 


Q's against the Z’s of the parents and as positive, and 8*+-decay Q's against the 
Z's of daughters and as negative. See Fig. 1(b). Alternate points are then 





Fic. 1. Illustrating plots of Q. The total decay energies are plotted against the smaller Z 
in the transition. Fig. 1(a) is for an odd A (135) and Fig. 1(b) for an even A (106). 


connected by straight lines. Two straight lines are obtained, one for odd Z 
and one for even Z, provided the spin-dependent terms refer to the same shells 
for all transitions. The slope of these straight lines is B, (equation [1.3]). 
However, in general, if any two points two Z units apart are joined, the slope 
of the line obtained will be 


[1.4] B’, = Bat 30 


where @ is the difference of the spin terms of the two Q’s. 
Moreover, we see from [1.2] that, provided no shell is crossed in transitions 


*The term will be used to include electron capture. 
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from one isobar to the next, the vertical displacement between the odd-Z and 
even-Z lines is 2(4; + »,) for even A and 2(7, — »,) for odd A. In these cases, 
we define 


[1.5] ba = wet VE and 64 * Fy — Vip 


Note also from [1.2] and [1.3] that if there are no shell crossings 
[1.6] Za = 3(Zi1+ Z2) +3 - (Py — m)/Ba — €4/Ba. 


Here Z, and Z2 are the intersections of the two lines with the Z-axis. The term 
e4/B, never exceeds 0.3 and is usually less than 0.1 (see §4). The term in p, — m 
may be of the same order of magnitude. Then the quantity Z’,4 obtained by 
adding 0.5 to the middle point of the segment of Z-axis between the odd-Z 
and even-Z lines will be a first approximation to Z,4. Fig. 1(b) represents one 
of the many examples studied. 

When there is only one line for both odd and even 2’s, i.e. e4 = 0, 2’, is 
obtained by adding 0.5 to the Z value at the point of intersection. An example 
is given in Fig. 1(a). 

The best available values of the various §-decay energies were taken, and 
a graph was plotted for each A for which there were data. Where actual Q’s or 
information on By coincidences were not available, considerations of shell 
structure, energy match, and consistency with neighboring decays were used. 
In a few cases it was possible to use Q’s from nuclear reaction thresholds. In 
still fewer cases mass spectroscopic differences, which are accurate enough in 
the low mass region, were used as auxiliary information. In the high-mass 
region Stern’s compilation of masses was used throughout (16). The survey 
of data (2, 10, 11, 14, 16, 18) was completed in the first quarter of 1954. Some 
points are liable to change as more reliable data become available. 

Wherever errors in experimental values were available corresponding limits 
were placed on By. Useful half Mev. wide limits can be obtained in cases 
involving electron capture. 

The effect of neutron magic numbers is such that the actual points corre- 
sponding to transitions involving these are found to lie above the liquid drop 
line, at Z corresponding to magic-plus-one neutrons. The proton magic numbers 
cause the corresponding point to lie below the liquid drop line at Z = magic. 
Such points are not generally used to determine B, and Z, values. However, 
at certain places, upper or lower limits can be placed with their help. 


3. Behavior of Bu 

In Fig. 2 are presented B’, values found by this treatment. It is not possible 
to determine B, or the spin terms quantitatively from the Q’s. (See equation 
[1.4]). However, using qualitative arguments about shell-binding and pairing, 
we can find whether B, is less than, equal to, or greater than B’,. Further, 
the same type of reasoning tells us for which mass numbers the values of By 
and B’, are closest. In Fig. 2 we plot those B’, values which are very near to 
actual B, values. The actual points obtained are put in circles; those with an 
arrow pointing downwards are upper limits. The limits obtained from the 
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Fic. 2. The parabolic constant By. The points are those computed from energy differences 
(§3). When limits are shown there is usually an electron capture involved of unknown energy. 


Curve a is fitted to these points, curve ) is obtained from packing fractions with variable a; (§8), 
and curve c is B, from the usual semiempirical formula with Fermi’s numerical values. 


electron capture data are also shown. The curve a passing through the points 
is the optimum curve and we consider it to be fairly close to the actual B,. 
It has been extrapolated at crossings of shells and in some other regions. The 
variations that appear on this curve may be due to spin terms or due to local 
variations in the average “‘liquid drop”’ function. 


4. Spin Terms 

It is not possible, in the present state of the data, to obtain these shell- 
dependent corrections empirically. However, we assume that e4 and 6,4 are 
constant for a given neutron and proton shell combination and that B, follows 
a smooth curve, such as a, Fig. 2. It is then possible to evaluate 6,4 and e, for 
various subshell combinations from the spacing between lines as in Fig. 1(d) 
and equation [1.5]. The transitions across the shells and on either side of magic 
numbers cannot be represented on this simple picture. For those cases we need 
to work with equation [1.1] and a fairly satisfactory correlation can be ob- 
tained, as will be seen in §7. 

The e4 and 6,4 values obtained are presented in Table I. The points from 
which these values were obtained are shown in Fig. 3. 

It will be noted that the crossing of mixed shells (e.g., g7/2, d3/2, $1/2, @s/2 for 
neutrons in region A = 106 to 138) does not affect the 6, term appreciably. 
It is reasonable to assume that after the sudden fall in region 138-156 6, will 
again become constant. A slight rise may be expected when 73,2 is crossed but 
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TABLE I 





Protons Neutrons 

3d5/2 4f7/2 

Afrye 4fije 

4fz/2 3p3/2 

4fzje 4fs/2 

3p3/2 4fsj2, Spij2 
3p3/2, 4fs/2 4fsj2, 3pije 
3p3p2, 4fs/2 5g9/2, 3pi/2 


3p3/2, 4fs/2 
3pi/2, 5g9/2 
5209/2 
5g9/2 
5g7/2 
5g7/2 4d5 2 


6hoj2, Sf7/2 


5g9/2 
4ds/2 
5g7/2 
6h41;2, 3512 
6hy {29 35172 
6 11/2) 351/2 
6hop, 5fz/2 
Ghop, Sfspe 
Tiisj2, 43/2 


\and 41/2, etc. 


54, Mev. €A, Mev. 


0.2 
0.0 
0.5 
0.3 
-0.3 
0.0 
0.0 
0.2 
0.5 
0.0 
0.0 

—0.3* 
0.5 
0.5 

t 0.0t 

0.3t 


mim mm NNNNNWNNN WN 
AnnnwwwwasSecooounts 





*e = Oafter A = 131. 


{Little evidence. 


SPIN TERMS Mev. 


150 


Fic. 3. Pairing energy sums and differences obtained from decay energies. 





is unlikely because of mixing with p1;2. Beyond A = 216 there is apparently 
a gradual decrease in 6,4. A constant value is, however, possible in view of the 
experimental uncertainty in the masses used (16). 

The situation is not very clear with respect to e,. Being a small quantity it 
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is very sensitive to any error in energies involved. Although the e,4’s are 
definitely non-zero, the magnitudes given here are not precise. Data are very 
rare in the 150-210 region. From analogy with 64 we have considered €, to be 
constant in the two regions 108-136 and 156-210. The terms are negative after 
crossing of proton shells and seem to confirm the rule given by Suess and 
Jensen (17).* These values may have to be changed somewhat as more data 


become available. 
5. Behavior of Z4 

In presenting information on charge numbers, it is convenient to use the 
quantity 7; = 4A — Z, which ishalf the neutron excess. It is seen from equation 
[1.1] that the minimum value of M(A,Z) occurs for Z = Z,4. We see from 
equation [1.6] that Z, = Z’, only if €4 and p, — m; are both equal to zero. 
Since the information on p, — m, and e, is not very abundant, we have been 
forced to deal with Z’, rather than Z,4. These values are converted to Tymin 
= $A — 2’, and this quantity is plotted in Fig. 4. Trm:, should lie within 
+0.5 of T; for the B-stable isobar. The values 3A — Zp for the stable odd-A 
nuclides are represented as horizontal straight lines in Fig. 4. The function 
Trmin is seen to be linear in A between shell regions. Changes of slopes and 
discontinuities occur at crossings of shells. In spite of different methods of 





a 
N 
1 — 
$ / 
= ? 
fe 49 
4. 
2 
- : 
a 
' 
: | 
et 
0 90 "Wo 130 150 OD 90 210 0 


mass-number A 
.FiG. 4. Tymin ='($A — Z,4) obtained from decay energies assuming the last two terms in 
equation (1.2) to be zero. The interpolations are made from the consideration of 8-stable odd-A 
isotopes. See §5 for discussion. 
*Each time a nucleon number crosses the magic number 28, 50, or 82, the pairing energy of 
this type of nucleon is smaller than that of the other type. 
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approach, these results are in general agreement with those of Coryell (4). It 
is not clear how his values will change in the light of our findings on B, which 
are generally smaller than the values assumed by him. 

It appears, however, that if it were possible to work with equation [1.6] 
and exact Z,’s the discontinuities across the shells might be removed and that 
Trmin Might yet be a smooth function of A. 


6. Classification of Information 

Apart from spin terms, we expect, on the basis of equations [1.2], to obtain, 
in general, two parallel lines on our plots such as Fig. 1. The lines on our plots 
are equivalent to isobaric parabolae. Since spin terms are different for different 
transitions in an isobaric sequence only when subshells are being crossed, we 
expect all the points to lie on the straight lines and parabolic dependence to 
hold when no subshells are being crossed. 

Evidence in this regard can be classified in three categories. To define a 
parabola we need three masses (two Q’s); to define two parallel parabolae we 
need four masses (three Q’s). If in such cases the B, values are not inconsistent 
with those from neighboring A’s, we can take it to mean that the parabola 
is ‘‘defined”’. We shall call such cases category 2. If, also, one or more addi- 
tional Q, lie on the lines so defined, we say the parabola is ‘“‘checked”’ and that 
these cases fall into category 1. If, however, there are Q’s lying away from the 
lines, the Z-dependence is not parabolic. This is category 3. In category 3 will 
be the cases where subshells are crossed. 

The cases in each classification are specified in Table I]. The cases where 
one line only is defined are those with e, ~ 0. There are 30 cases where 
parabolae are checked, 16 where they are just “‘defined”, and 15 where the 


TABLE II 
CLASSIFICATION INTO CATEGORIES 1, 2, 3 


— 








Straight lines 


No. of cases No. of Q’s defined A 
Category 1 
2 3 1 109, 115 
2 3 +EC 1 111, 113 
2 4+EC 1 135, 133 
11 3+EC 2 {73 82, 83, 85, 101, 105, 
127, 129, 134, 151, 192 
4 4 2 96, 99, 106, 112 
6 3* 1 55, 57, 71, 73, 75, 79 
3 2+ 2EC’s 1 97, 171, 161 
Category 2 
4 2 1 54f, 170, 166, 150 
12 3 2 (i 56, 70, 76, 78, 86, 94, 
130, 132, 144, 146, 156 
Category 3 
2 4*t 1 107, 77 
6 4 2 49, 84, 87, 89, 92 
2 3 + EC*t 2 146, 144 
5 3 + ECt 2 88, 93, 137, 139, 140 


*Small spin effects present. 
tAssignment doubtful. 
} Magic numbers involved. 
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Q2+2+, Mev. 





' 
a] 


-4 


-6 


x estimated 
(Sec.7) 
aw uncertain 


Qz+z4, Mev. 





Fic. 5. Specific cases illustrating use of equations [1.1], [1.2], and [1.3] when shells are 
crossed. See §7. 


parabolae are not followed. These last 15 cases occur mostly where major 
shells are being crossed and in the rest of the cases also, spin changes across 


subshells are occurring. It may, then, be said that the parabolic dependence 
is correct, except as modified by spin terms. 


7. Specific Cases 


We have seen that the parabolic law breaks down when shells are being 
crossed. A simple explanation of this effect can be given on the basis of equations 
[1.1] and [1.2] and the shell structure of the nuclei involved. According to the 
shell model we expect the pairing energies to increase with /-value and for the 
same /-value the larger j-value should have the greater binding energy. 
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The cases we discuss here are those at the crossing of the 50 neutron shell. 
Cases occurring at the crossing of other shells show a similar pattern in their 
energies. We use the level scheme given by Klinkenberg (13). The neutron 
shells involved in the transitions are 4g9/2, 4ds5;2, and 4g7/2. We expect v9/2 
= 7/2 = V5/25 and N9/2 a N5 /2 2 N7 /2- 

The proton shells involved are 3f5/2, 33/2, 31/2, and 4g9;2. We expect 
19/2 > 5/2 > 32 = TMiy2 and ps5;2 = p3j2 > piy2 = Poo. On the basis of 
these relations and rough estimates of the order of magnitude involved, we 
can compare the various spin terms for nuclei 87 < A < 91. The data on 
which this discussion is based are presented in Fig. 5. 


A = 87 


Q35 is evaluated from equation [1.1], using the fact that, in the nucleus 
Z = 35, the 35th proton is in a 3/2 state, while in the daughter nucleus the 
35th and 36th protons form an fs,;2 pair. The remaining Q’s are found from 


[1.2] as: 
35 = LQ35+2)s)2— P32—52+2 45/2 — 25,2, 
Q36 = LO36+P3/2—N5/2, 

7 = LO3+P3/2—9)2+2 43 /2— 29/2, 

3 = LO3+P12—N9/2, 

39 = LQ3+P12—M92+2 41/2— 29/2. 


The B’, values from alternate points in Fig. 5(a) are 
I g 


B's1g = Baz + (bs/2 — P32) + 3(M9/;2 — M52) + (45/2 — 73/2) + (9/2 — ¥5/2) 


= 3.8 Mev. from Z = 35, 37, 
B’ sin = Bs; = 3[ (P32 — P1/2) + (N92 = N52) = 2.9 Mev. from 36, 38, 
B's, = Bgz + 4(p3/2 a P1/2) + (13/2 —_ 11/2) = Mev. from 37, 39. 


From a consideration of the relations between the various spin terms explained 
above, we expect B’s7, > B’s7, > B’s7- = Ba, which is qualitatively correct. 
The difference between cases a and b, however, seems rather large; we shall 
return to this after examining A = 88 and 89. Of course, if Bsz were known, 
some of the spin differences could have been estimated. This case illustrates 
the large effect of a change of the state of the last particles. 


A = 88 





In these cases, equations [1.2] give 

36 = LQse+P3/2—M52—2¥5/2, 
37 = LQ3+P32—Ms2+2 73/2, 
33 = LOQ38+P1/2—M9/2— 25/2, 
39 = LQ3+p12—M92+2 71:2. 


Therefore, the B’,’s in Fig. 5(6) are 
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B' seq = Beg + 3[(b3;2 — piy2) + (m9/;2 — M5/2)] + (9/2 — ¥5/2) 


= 3.4 Mev. from 36, 38, 
B’ ssp = Bas + 3[(b3;2 — pij2) + (9/2 — M5/2)] + (43/2 — 1/2) 
= 3.1+0.25 Mev. from 37, 39. 


As before, we expect B’ssg > B’ss,, which is verified. Moreover, since 
Bss = Bgz, inspection of spin terms leads us to expect B’s7, to be slightly 
greater than B’ss, and B’s;, = B's). The experimental values, from Fig. 5(a) 
and (6), are 3.8 Mev. and 3.4 Mev., 2.9 Mev. and 3.1 Mev., which fit in well 
with these predictions and with the orders of magnitude of the quantities 
involved. 


A = 89 
The Q’s from equation [1.2] are 
89 89 
Q3s = LQset+ps/2—M5/2, 
89 89 ; 
Q37 = LQ37+p3/2— 5/2 +2 13 /2— 25/2, 
89 89 
Q3s = LQ3s+fij2—M5/2, 
89 89 9 9 
Q39 = LQ39+P1;2—M9/2 +2 1 /2—209/2. 


Therefore, the B’,4’s (Fig. 5(c)) are 


B’ 394 = Bgg + 3 (p3/2 oe P1/2) ao Leg Mev. from 36, 38) 
B’ s9n = Bsy + 3[(P3/2 — P1y2) + (m9/2 — M5/2)] + (43/2 — miy2) + (¥9/2 — 5/2) 
= 3.9 Mev. from 37, 39. 


We expect that B’s9, > B’s9g. This again is verified. The value Bg, is very 
small on the basis of the present data. However the decay scheme of Kr°® is 
not well established. 

A further confirmation of this type of reasoning comes from the following. 
From the spin terms involved we expect B’s7, > B’s9, > B’ssq. The values 
actually obtained are 3.8, 3.6, and 3.4 Mev., respectively. 


A = 90 
The Q’s from equation [1.2] are 
Qis = LQ36+P3/2—M5/2— 25/2, 
37 = LQs+Ps2—M52+279/2, 
38 = LO3s+P1j2—Ms)2—2y5 2, 
39 = LOQ3t+Pip—Ms2+2m1)2, 
Qto = LQ0+Ps2—M92— 29/2. 
. The B’ ,’s (Fig. 5d) are 
B'soq = Boo + 3(3/2 — piyz) = 1.3 Mev. from 36, 38) 


B' sop = Bgo + 3(ps/2 2 P12) ae 73/2 — Tij2 = ee Mev. from 31, 39, 
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B’ ye = Boo + 3[(Piy2 — Pose) + (M9/2 — 5/2)) 
+ (v9/2 -— 5/2) = 2.6 Mev. from 38, 40. 


In spite of slight uncertainty in the data, these values are also in keeping 
with the arguments given before. 

A=9l 

On the basis of existing experimental values this case can not be fitted into 
this kind of argument. However, the energy measurements are by no means 
complete. See Fig. 5(e). In the decay of Kr®! and Rb®! only the 6’s have been 
reported; the y’s observed were not measured. We may ‘se the better known 
energies in the decays from cases just discussed to pre’! the energies of the 
doubtful decays. 

The Q's are, from equation [1.2], 


91 91 
36 = LQ36+P3/2—N5/2, 
91 91 9 9 
37 = LQ37+P3/2— 5/2 +2 13/2 — 25/2, 
91 91 
Q33 = LO3s+fijy2—M5/2, 
91 91 9 9 
Q39 LQ 39+ P12 — 5/2 +2 m1 2 — 25/2, 
91 91 
Qio = LQwtPo2—M5/2, 


91 91 ’ 
41 LQa+P9/2— 19/2 +2 19/2 — 29/2. 


The B’,’s are 


B’sig = Bor + 3(b3/2 — Pipe) from 36, 38, 
B’ yin = Bor + 3 (P32 — pry2) + (43/2 — 1/2) from 37, 39, 
B’ sic = Bor + 3(Piy2 — Poy) + (1/2 — 79/2) 

+ 3(n9/;2 — M5;2) + (v9/2 — ¥5/2) from 39, 41, 
B' sia = Boi + 3(Pis2 — Po/2) from 38, 40. 


Assuming Boo ~ Bog; we have 
[1.7a] B' sig = B’ soa, 
[1.75] B’ sip = B’ soo, 
[1.7c] B’ sie < B’ soc. 


From [1.7a], the well-known decay of 3sSr®*', and values from Fig. 5(d), we 
find that energy in the decay of 3sKr®! should be 5.3 Mev. The value is at 
present reported as >3.6 Mev. Similarly, from [1.7b] and 39Y°"’s well-known 
decay, 39Rb®! should have Q = 5 Mev., which is to be compared with the 
experimental value (> 4.68). However the comparison of B’s;, and B's, 
suggests that 9/2 — 1/2 is almost zero, which is surprising. 

In this manner, with the help of known energies and this analysis, a satis- 
factory estimation of unknown energies and, therefore, masses across the shell 
regions can be made. A few cases may be mentioned. 5:Sb'** and 52:Te'*® are 
known to be 8--active. A determination of these energies would be helpful. 
From Fig. 1, the predicted Q values are 4.1 and 3.1 Mev., respectively. 5sCe'** 
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is expected to undergo 6+-decay with a Q of 3.4 Mev. Experimentally, a 
8-spectrum with 1.3 Mev. anda y-ray of 1.8 Mev. have been found, with no 
information about coincidences. Sb'*° are known to be 8--active. Our extra- 
polated Q’s are 1.1 and 4.2 Mev. For studying spin terms near shells, the Q’s 
of the following B--decays would be valuable: (I, Xe, Cs)'89 and Cs!8°, 


II. ABSOLUTE VALUES OF MASSES 


8. The Liquid Drop Constants 
Here we try to fit the information obtained in Section I within the frame- 
work of a liquid drop formula. The usual liquid drop formula is written 
[2.1] M(Z,A) = M,A — (My — My)Z — iA + a2A*? + 43(3A — Z)?/A 
+ Z*cA-'/3 + spin terms 


where V/, is the mass of neutron, My is the mass of hydrogen atom, a@;A is the 
volume energy, a@2A*/* is the surface energy, a3 is the symmetry energy con- 
stant, and c = 3 e?/ro is the coulomb energy constant. The nuclear radius is 
roA'/3 and we define A = M, — Mg. The spin term is assumed to be zero for 
odd-A and +36, for even-A, odd-Z and —36, for even-A, even-Z nuclei. This 
approximation is justified in view of the finding that €4 is quite small for most 
values of A. 

From Equation [2.1] we get the charge of the most stable isobar by putting 


ies oy = 0, giving Z, = a: 
“(03 

If we put 

2.3] By = 2(a3 + cA?”)/A 


we can write 1/(Z,A) = M(Z4,A)+ 4B4(Z — Z4)? + spin terms. 
We have obtained Z,4 and B’, values in Part J. If the liquid drop formu- 
lation is to be valid, then from equations [2.2] and [2.3] we should have 


[2.4] Z4B4 = a3 + A = constant 


as has usually been assumed in earlier works. Even taking into account the 
fact that B’, and By, are not precisely equal (equation [1.4]) the variations in 
the product are not compatible with a constant a3. Equation [2.4] is inde- 
pendent of whatever values the other constants may have and it thus forces 
us to conclude that a; is not a constant throughout the mass table but varies 
from mass-number to mass-number. 

Alternatively, we may eliminate a; from equations [2.2] and [2.3], obtaining 
[2.5] c = 4B,A'* — (Z,B, — A)A-*'. 


In this case we find, using experimental values, that c is not constant. Since 
this factor is used to define the coulomb radius we can examine what effect 
its variation would have on ro. Using the experimental By, values of Fig. 2 and 
the natural limits on Z,, i.e. Zp+4, we have calculated the limits on c from 
[2.5] and present these as values of ro in Table III. These large random vari- 
ations in ro are in direct contradiction with the evidence from many other 
sources that 7o is quite constant (e.g. References 1, 5, and 15). We conclude 
that variations in the factor c cannot be used to account for the discrepancy. 
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TABLE III 


LIMITS ON COULOMB RADIUS CONSTANT fo 
CALCULATED USING By, OF F1G. 2 AND NATURAL LIMITS ON Z, 


ro, X 1078 cm. 








| 6. 6 10-8 cm. 





| 
| 
A Upper Lower | | A Upper Lower 
53 2.37 1.81 || 135 1.99 1.84 
57 2.57 Ls || 147 2.17 2.00 
71 1.37 112 |) 151 1.65 1.53 
73 1.61 131 || 171 1.4 1.31 
77 2.19 1.79 || 185 1.22 1.15 
81 1.86 1.57 || 213 1.02 0.97 
105 1.76 1.52 || 223 1.15 1.1 
111 1.86 1.64 || 225 1.31 1.25 
113 1.53 146 || 231 1.4 1.346 
237 1.24 12 


The alternative explanation is that [2.4] and [2.5] are incorrect, owing to 
some invalid simplification in the Z- or A-dependence of [2.1]. We should 
point out that this may be expected since we are here trying to fit individual 
masses accurately rather than to obtain an agreement between [2.1] and an 
average trend of nuclear masses. 

As a simple generalization of [2.1] it may be possible to keep its form but 
to let just one of the ‘‘constants’’ vary with A. To explore this possibility, we 
have attempted to put all the variation in a3. For the packing fraction of 
stable, odd-A elements, the odd-even terms are zero by definition and the 
parabolic-terms are quite small. Hence, for the packing fraction (14 — A)/A 
we put M = M(A,Z,) and obtain 


[2.6] P.F. = -a+a,A7?%+ 

where —a = 4(M, + Mag) — (a: + 1000) m.M.U. 
and = = ic(a3 + 2A)A?/3(a3 + cA?/%)-1, 

We note 

[2.7] a3 = cA*/3(4E — 2A)/(cA?/3 — 42). 


Making a suitable choice of the constants c, a, and a2 and using experimental 
packing fractions of odd-A 8-stable elements, we calculate a; using equations 
[2.6] and [2.7]. 

Since we have chosen a; for each A to give the experimental packing fraction, 
the numbers which remain to be fitted are B, and the natural limits on Z,, 
viz. Z> +4 > Z4 > Zo — 3. Of course, ideally, Z4 should agree with those 
values obtained in Section 5. The variations of a3 (and therefore of B, and Z, 
thus calculated) depend for their magnitude on the choice of other constants. 
The manner of variation is, however, little affected for any set of constants 
which give a reasonable fit to the data. The values of B, obtained in this way 
from the packing fractions are shown in curve b of Fig. 2. It is seen that the 
agreement with the values deduced in Part I from isobaric mass differences is 
not good. Also for comparison we show in curve c the values of B, from the 
liquid drop formula [2.3] using the constants given by Fermi (quoted in 
Reference 8). 
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The value of c can be fixed within rather narrow limits. In earlier works 
values of c from 0.62 to 0.75 m.M.U. have been used. The points in favor of 
using a high value of c are (1) closer agreement with recent values of ro (1, 3, 9, 
12) and (2) less fluctuation in Z,4 (calculated as above). The points against a 
high value are (1) generally smaller values of Z4 compared to natural limits 
and (2) generally larger values of B4 compared to curve a (Fig. 2). 

Atc = 0.75 m.M.U. calculated Z,4 values were 3 to 4 charge units too small. 
For higher values the difference increased still more. At c = 0.7 m.M.U. the 
deviations were reduced to about 0 to 1.5. In Fig. 6 we show Z4 = Zo — Za 





Fic. 6. Zo — Z4 obtained from packing fractions. The points should lie within the dotted 
lines. Values are shown for c = 0.686 m.M.U.(0) and for c = 0.7 m.M.U.(@) for nearly the 
same values of Bg. 


for the two choices c = 0.700 m.M.U. and c = 0.686 m.M.U. It is seen that 
the latter value is preferable. We cannot go below 0.686 m.M.U. because the 
fluctuations in Z,4 increase and with them the number of disagreements. Other 
values of a and a2 were tried without improvement. For the points corre- 
sponding to c = 0.686 m.M.U. in Fig. 6 the values are a2 = 16.33 m.M.U., 
a = 7.3526 m.M.U. (or a; = 15.91 m.M.U.). 

These values are, however, no contradiction of earlier works. The disparity 
arises because we are taking individual experimental masses whereas in other 
works a departure in masses to the extent of 5 m.M.U. was tolerated. It is 
interesting to note that a departure in masses of the order of 3 m.M.U. is 
sufficient to bring both Z, and B, to reasonable values for medium values of A. 
As A increases, we need still larger departures from actual masses to remove the 
stability shifts. 

The fluctuations in stability shift (2) — Z4) shown in Fig. 6 are alternately 
in opposite directions. This makes it impossible to improve the broad range fit 
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by any adjustment of constants. While the agreement in any two adjacent 
sections (Fig. 5) can be improved by a suitable choice of constants, it becomes 
worse in the contiguous sections on both sides. 

To correct for the local stability shift, we observe that at places where Z, 
is about a unit away we need about 2.5-3% change in ¢ to bring Z,4 to the 
nearest natural limit. The same thing can be achieved by a 2% change in a2 
or a 0.8% change in a. No attempt is made, however, to go any further in this 
direction because of the arbitrariness involved in such a procedure. 

We are thus forced to conclude that varying just one constant a; is not 
sufficient to account simultaneously for the various fine structures in both the 
Z, and packing fraction. It appears that the A- or Z-dependence of the formula 
is more complicated than the one usually assumed. In this connection, we are 
examining the effects of non-spherical nuclear shapes which affect the order in 
which single-particle levels are filled and may cause important corrections to 


the masses. 


9. Empirical Computation of Masses 

The first requirement in computing masses is to know the mass A/(A,Z,4) 
of the hypothetical nucleus with charge Z,. We have seen in this work that 
this cannot be obtained very accurately (say, within 1 m.M.U.) without 
introducing some complicated A-dependence in more than one “‘constant”’ of 
the mass formula. For purposes of calculating masses, we present an empirical 
interpolation method. 

We know that masses of stable odd-A elements are very nearly equal to 
M(A,Z.4). We have therefore plotted in Fig. 7 the odd-A packing fraction 


odd A 
evenA 
Gorregted 
even 


(Pfexp— Pfr) * 10° for “stable’ atom 





Fic. 7. The packing fraction differences AP.F. defined in §9. Even-A points are corrected 
for 54. 
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differences AP.F.(=P.F..xp — P.F.tn) against A. P.F.,, could be any con- 
venient reference level and the values we have chosen are those calculated 
from a liquid drop formula where all constants (including a3) are fixed at the 
following suitable values: 
a; = 16.89 m.M.U. (a = 8.3341), a2 = 19.45 m.M.U., a3 = 100 m.M.U., 
c = 0.75 m.M.U., A = 0.84 m.M.U., 


P.F..n = —8.3341+19.45A-1/3+-0.19066A ?/*/(1+0.0075A?/*) for odd A’s. 


In regions where odd-A masses are not known, even-A AP.F.’s are also 
plotted. A comparison with Fig. 2 reveals that variations in AP.F.’s show 
even more marked regularities with A. The large dips occur at the crossing 
of major shells. This information is useful since it helps fix the trend of the 
curve from analogy with known parts. 

In the region 135-200 only a few odd-A masses are known and the curve has 
been fixed using even-A packing fractions. To correct even-A masses we use 
the 6, values of Table I and the approximation involved therein. The approxi- 
mation is good enough in view of experimental errors involved in masses in 
this region (shown in Fig. 7). 

We have tried to draw a smooth packing fraction curve. The masses taken 
from it may, tentatively, be considered good within 1 m.M.U. or less. The 
amount of experimental error is shown in the high mass region; in the low mass 
region it is contained within the circles. The whole family of Stern’s masses 
203-245 can be shifted either way by the amount indicated on the graph (16). 
Although nothing can be said very definitely, it appears that an upward shift 
may not be out of place. 

For the purpose of computing masses in regions where shells are not crossed, 
equations [1.2] and [1.3] can be used with values of Z, taken from 7;m1, curve 
of Fig. 4 and ignoring p; — n;. The dotted lines show the tentative interpo- 
lation where 7; m1, cannot be obtained from energies. B,’s should be taken from 
curve a in Fig. 2 and spin terms from Table II. 

However, where major shells are being crossed the simple approximation is 
not very close. In such cases equation [1.1] can be used with a fair amount of 
success as the spin terms involved frequently give similar differences, and 
information from a well-known A can be used in the neighboring A’s. We have 
illustrated the method involved by giving a few specific cases in Section 7. 
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LONG DISTANCE V.H.F. FIELDS 
II. REFRACTIVITY PROFILES CONTAINING “SHARP LAYERS”! 


By F. H. NortTHOVER? 


ABSTRACT 


In this part the general solution developed in Part I is applied to the case in 
which the refractivity profile contains sharp layers (i.e. local regions where bu 
and dyu/dh change very rapidly compared with their rate of variation in a “‘stan- 
dard”’ type atmosphere). It is found that such layers, when well developed, can 
cause distant fields of the order of magnitude of those which have been observed, 
but present experimental evidence seems to indicate that scattering from atmo- 
spheric turbulence is usually the important factor. An attempt is made to work 
out a physical interpretation of the field-formation below the elevated layer. 


SECTION 1 
1.1 Case of One Layer 
We shall suppose that, in a typical profile, du/dh changes only slowly with 
height nearly everywhere, the exceptional regions being localized in narrow 
bands where this change is relatively very rapid. These bands we shall refer to 
s ‘sharp layers’, or, briefly, as ‘‘layers’’. As has already been remarked, they 
occur quite commonly at heights in excess of 5000 ft., especially in the height 
range 5000-10,000 ft. We shall consider first the case when only one layer is 
present and then try to see how the theory can be extended to take account 
of the case when multiple layers occur. 


When 
(63]" Aamitle—et Bs 


is small (case of profile of §4, Part I), then 





3Q° 19" 5g , 9m" 
[64] 40° 20 36W ~ 7 (—m) ats (i—m)'s ’ 
while if 

: m' 
[65] ina 36| > 1, 

$C if so 3 elec 

[66] 40°20 RW ip 
Integrating [41] and putting p = 0, we have 

7 er ses ge 
[67] Fo) = Jr 1020 36W/)F’? 





and only from local regions where du/dh is changing so rapidly that condition 
[65] is fulfilled is contribution to this integral likely to be sensible. 


1Manuscript received February 1, 1955. 

The work reported in the present paper was carried out in the laboratories of the National 
Research Counci! in consultation with N.R.C. staff during the summer months of 1954. 

Department of Mathematics, Memorial University, St. John's, Newfoundland. 

*This numbering follows on from that of Part I (Can. J. Phys. 33: 241. 1955). 
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1.2. Case of One Elevated Layer 

If the layer is sufficiently thin* we may, since Q(p) is known to be a continu- 
ous function,f take its value throughout the layer as being practically equal 
to its value at the upper boundary. Let the height of the layer be Rp;; then, 
throughout the layer we may take, effectively, 
[68] Q2(p) = F(p1) 


since 
@ 3 72 r 2 1 ” Q 
a) = FW)- ['EG-3 35) pa 


and for p > p; the integral is of order m,?/Z,f (the subscript 1 denoting values 
at p = pi). 
Hence, throughout the layer, 


[69] Q(p) = ilexp{ —$V/2.ix(p1— e+ 361)*? 1} +7] 

if 71 x7/3(p;—e+46;) > 1, mw > arg(pi—e+45;) > 4x 
(here y; = 1/(1—m)); 

[70] Q(p) = iexp{ —44/2.ix(p1— e+ 361)*? y1} 

if x?/3(p;—e+46:)y1 > 1, te > arg(pi1—e+36;) > 0; 


and Q(p) equals 
Hy iV 2.x (pret $1)" 91) 





[71] HV (2/2.x(o1—e+ 301) vs} 

if x?/3|o,—e+36;| < 1. 

Also 

[72] F'(p1) ~ 2/2. (p1— €+ 361)? exp|{ —$V/2.ix (p1— e+ 361)*? v1} 

if x*/3(p;—e+461) > 1 

and (by Appendix 3) 

73] Fp) = —Se**—— aT 445,751] 
” 7 yi (pi—e+ $51) [AY (3/2. (p1— + 451) 1} 

if x?/8|o,—e+36:| < 1. 


Hence, if (Am), denotes the total change in ‘‘m”’ as we proceed through the 
layer, the approximate value of the contribution of the layer to the integral 


on the right of [67] is 
= +e (Am), — ° =e ‘ ‘a ra 1 3/2 
[74] 8/2.x(p1—e+ 90s) 7 exp{ —3V/2.1(p1—e+ 351) 1} 


x [1 +7 exp {$+/2.xi (o1—e+ $61) "'"71} ] 
if x?/3(p,—e+46:)71 > 1, mw > arg(pi—e+46;) > 4x; 


*A ppendix 6. 
tSee §4. 
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=e (Am); 4 3 aes 
sigs = a ca 15 

(75] 8./2.x(pi—e+-3 15, exp{ —3V/2.x1(91-e+ 361) 1} 
if x*/3(py—e+ 361) 71 > 1, $m > arg(pi—€+361) > 0; 
and 
[76] — 37 (Am), eH? (5/2. (o1—e+ $51) rs} 
if x?/3|p,—e+46)| Se i. 

It is easy to verify that the asymptotic expansion of [76] gives [74] as a first 
approximation. 


Since contributions to the integral on the right of [67] outside the sharp 
layer are negligible (being at most of order mo?/Zo) it follows from [66] and 
[67] that the approximate equation for e is 


HY (5/2. x(—e+4 159) * 7-0} 








7] PV 2x(— Be) yo) 
A : 3/2 
(a 15,)° exp {3i7— 3V2.ix(p1—e+ 461) v1} 
when x?/3|o,—e+46,| > 1, 0 < arg(pi—e +361) < 3x 


which is exactly the same as the equation obtained by Carroll for the bilinear 
model (Reference 3, p. 16, equation 24).* This is as it should be, for any linear 
portion of the profile will give a strictly zero contribution to the height integral 
on the right of [67]. We thus obtain the bilinear model as a special case of this 
theory. 

In practice, the great difference between the sort of sharp inversion layer 
physically possible and Carroll’s model lies in the fact that the former can 
occur within quite a wide height range, whereas the height of Carroll’s du/dh 
discontinuity is fixed and near the probable upper limit of the heights at 
which the inversions can in practice occur. 


Note.—A more exact form of [77] is 


7a) HPV 2(—et 40)? rol 
HOE SV 2.x(—e+ bo) yo} 


= —!2%1(4m) (1H [1 ? GV 2.x (p1— e+ 361) 7971} ° 
which holds when 
x?/3|o,—e+36:| <1. 


It will, however, transpire that roots for which [77A] would have to be used 
give negligible contribution to the field strength. 


2. EVALUATION OF THE FIELD FUNCTION GIVEN BY THE 
EIGENVALUES 


By equation [10] the field function consists of the sum of terms of the form 
F ,0, where 
[78] F wd Uy—4(b ) uy (E) 
J 





| uy -4 (x) { Ous- 3 (x) / OS} (sm) : 
*His s§ = 1—(Am), and his Am*/? = xv/2. (p1—e+48))9? 41. 
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[79] 0; = v,;P.,-,(—cos 8) . sec vz. 


The factor F,; depends on the heights of transmitter, receiver, and inversion 
layer but not on angular distance 6, while 0, gives the attenuation of the 
various modes with distance R@. We shall refer to the factor F, as the ‘‘ampli- 
tude”’ factor. 


Evaluation of the Factor F, 
We have (the dash now denoting d/dé) 
uy_y(x) = (Ay~itBy2)’ 


= (A V22) (pmo) 
= A(0)¥2(0)F’(0) 
27A (0) 
80 = ——— 
[30] ¥:(0) 
by Appendix 3. 
Now, for p > 1, 

[81] Q(p) = F(p)+k 
where 
[S14] k= Spam): eT TAD {3/ 2.x (p1—e + 481) rat. 


From the form of F(p) it is clear that 


OF 1 oF 
9 A a gt 
[82] Os l—m 0&° 
Hence 
: ay 2 ee *) 
[83] ( Os ) ~ l—mo 2h ( Os J 6 


(using Appendix 3) 
sah gee s a f ( ax} 
" (l-m) woe UE 21(1—amo){ ¥2(0)} ds J )* 
Noting now that 
[84] { Ou s_4(x)/AS} (ser) = A(O)Y2(0) (AQ/ds)o 


we have 


Uy—4 (x) { Ou s_4 (x) / AS} (s—) 
; ~. — 41401" _f,_ 3, ( *) \ 
[85] a (i—m) {20772 3t(1 mo) { ¥2(0)} as ye 
To evaluate the factor u,_,(b) u,_4(&) note that, below the inversion layer, 
du/dh is changing relatively slowly and consequently (see §4, Part I) A’ is 
small of order mo?/Zp, i.e., |A’| K 1. Hence we may, below the inversion layer, 
take A as effectively constant, so that, approximately,* 


* Assuming that & ts the affix of a point below the inversion layer. 
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[86] u,-4(b)u,_4(—)= {A (0)}? Q(pr)Q(pr) Y2(or) ¥2(pr) 

where pz, pe are the values of p at the transmitter and receiver respectively. 
Hence 

[87] F, ay 1d —m is (or) +kYs(or) | {ya (or) +RV2(oR)h yy (gy)? 


—4i(1—mo){ ¥2(0)}"(dR/ds)o 


{vilpr) ¥2(0) — ¥2 (or) ¥i(0)} { ¥i(pr) ¥2(0) a= ¥2(pr) ¥1(0)} 


= (1—mo) {1—41(1—mo) { ¥2(0)}?(d2/as) 0} 


and, of course, 


0 ix 2 < ‘ 2 
S)- Ty (Amn) eG? iy/ 2.2 (pet $61) °? 40} V2. (01 e+ 961) 


OS i 
‘ 7 (2) 2. a 15, 3/2 
x [ eR Rented) My) Evasion | 
88) ~ ~ EG nck fa) OPI BV 2a ret) 


since, in the applications with which we are concerned, 
R \x?!?(py—e+461)} D1 
with SF (xe) = O(1). 

For purposes of comparison with the “effective radius” solution, note that 
the latter is regained if, throughout the preceding analysis, we make (Am), 
equal to zero, for then k becomes identically zero. 

The field is now determined once we know the eigenvalues (the permitted 
values of «). These are given by equation [77]* and must be carefully studied 
before further progress can be made. This study will be the subject of the 
section which follows. 


3. THE EVALUATION OF THE EIGENVALUES 
3.1. Preliminary Consideration of a ‘‘Typical’’ Layer 
Fig. 1 shows the sort of profile which we might, perhaps, take as reasonably 
representative of an elevated inversion layer at, say, 6000 ft. 
> 
‘ 
& 


ide de nee sem ene 





Height 


° 
© 
uw 
a 
2 


*Or by [77A] if appropriate. 
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Our mathematical analysis has shown that the significant portion of the 
layer is that part of it where m’ is large (the ‘‘onset’’ AB of the inversion), and 
this is a relatively small portion of the “‘observable’’ layer* AC. 

If we take AC to be 100 meters thick (this is not impossibly thin) it would, 
perhaps, be reasonable to take AB as 20 meters thick. In that case, 54—ds 
would probably be about }(54—6c). Now in the “‘Cardington” inversion (10), 
64—5¢ = 5.16X10-* and this is quite large (although there have been high 
level inversions with 64—4¢ as large as 8X10-5). Taking Carroll's value of 6o 
as standard, it would probably not be unreasonable to take 64—5c¢ = }; do, 
which gives 64—5p = 4 59. Then the average value of m’ through the onset 
AB of the inversion is 


~ oF —4 . 9\2 
1 7.26 X10 = 3.63 X (3.2) «x 10°. 


2x me 20 F 6 
6.4 x 10° 


Now (see Sections 3.3 and 3.4) the eigenvalues which give any significant 
contribution to the field make 
ler — e+ 363| > x78, 


Even for \ as small as 10 cm., we have x~?/8 = (1/1287)?/3 10-4, so that all the 
significant eigenvalues make 


, 
a e myiler— e+ Bah > 1 
and the condition [65] is satisfied. 

The inversion layer can be regarded as a “local deformation”’ of the profile 
from the slow, steady decrease characteristic of the standard atmosphere type 
of profile. From this point of view, the values of m just below and just above 
the inversion must be taken as practically equal and equivalent to the value 
to be expected at this height with a ‘‘standard’”’ type atmosphere. We shall 
denote this value of m by m,, to correspond with the inversion height Rp. 
For example, for the exponential type of standard atmosphere, 


[89] mm, = mo exp (—2mopi/6»). 


Naturally, we write also 


[90] v1 = 1/(1—m)). 

Let 

[91] x?/3y79?/3 (e— 359) = X—iY, 
[92] x%Bq == Zo, 217/98, = Z). 
Then 

[93] x(p1—€+ 361)*? y1 


= [n?{ Pi} Zore**—Z1)} — (X-4Y) (1n/70)?? 


where P, is the value of P (see equation [44], Part I) at the inversion, at which 
p = pi. 
*I.e. the part where the profile slope is large. 
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Let 
[94] A = y173{Pi—3 (Zoro? —Z;)}, 
then we shall call A the ‘‘modified” or ‘‘effective’’ value* of P:. The quantity 
RAx-*/3 could then aptly be called the ‘‘effective height’’ of the inversion. 
Then, since mp and m, are both small compared with unity, we may write 
(y1/yYo)?/* = 1, and we have 
[95] x (pi—e+46))9 71 = (A—X+iY)3”? 
except when the right of [95] is < 1. This case, however, does not have to be 
considered for it turns out that eigenvalues for which the right of [95] is small 
give negligible contribution to the distant field. 
Hence the equation in for the eigenvalues is 
Hy” 54/2.(—X+iY) T= L (2) 3/2 
=— Am) [Hy (3./2.(A—X+iY 
H,” fz J/2.(—X+iV)*” * = al i3V ( +1 ) iy 
yi(Am)1 
“8 2-245)” 
For the particular layer we have been considering above, the average value of 
m within the layer is evidently about 2. We shall therefore take the value 
2 for (Am). 


[96] 5 5 exp{jim—jv/2. (A—-X+i1Y yay i. 


TABLE I 
REL ATION BETWEEN Pi AND A FOR io = 10 CM. OR A = 6 METERS 








Height P; P, Py A A A my ¥y1 
(ft.) (\=6 meters) (A=10cm.) (A=6 meters) (A=10 cm.) 

7,623 0.525 12.90 107.75 TAOP; 14.19 217-5 . 1947 1.24 
11,435 0.7520 19.35 296.63 1.029P, 19.90 305.1 Adi8 1.206 
15,246 1.00Z5 25.80 395.5 0.9936P; 26.63 392.9 o16.. L228 
19,058 1.2525 32.2) 494.4 O.9720P, 31.35 480.7 .13388 1.154 
22,870 1.50Z5 38.70 593.3 0.9593P, 37.26 569 .0 .1181 1.134 
26,682 1.75Z¢ 45.15 692.2 0.950P, 42.91 657 .7 1042 = 1.116 
30,4942 2.0025 51.60 791.1 O.944P; 48.71 746.7 .092 1.101 


“At th is height we reach the du/dh d iscontinuity mn Carroll's model. : 
3.2. The Equation for the Eigenvalues 

This is equation [96] above. As in the writer’s previous paper (10) (where a 
similar equation is encountered) the discussion of the equation falls naturally 
into the three cases X K —1, X = O(1), and X > 1. 

Case]. pk Sat 

Equation [96] becomes, approximately, 


* A is nas ai 
[97] -l= 8/2 TRS San exp{ — —3/2.(A—X +7) Ob 
és; 
et i nel ha ee ia 


8/2.(A—X+iY) 


provided X is not too near A, for if X is not too near A 


*In the case of an airless earth we should have A = P,. 











~ 
to 
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[98A ] Y/(A-X) «1 
owing to the large size of P; for inversion heights and wavelengths likely to be 
of importance in practice. Hence, approximately, 
—3V/2.i(A—X)*?42{2(A—X)}°"_¥ 


‘ 3/2 aia 
[99] = —(2p+1 Jin Hog SAS E : ret 
Equating imaginary parts 
[100] A—X = 30), approx. 
where 
o,/? = 3(2p+1)r 
[100A } = 3(p+}4)r. 
Equating real parts and using [100], we have 
3/2 
a Y= 3m oe( Ae) 


provided, of course, that 
l 4a," 

375 | + a 

732 og (ies, <K 
1.e., that 

1 12(p+3)x 
102 = log( 22 —2--} & 1, 
[102] 3(pth)e 8 y1(am); 
In the case we are considering, (Am), ~ 2, so that the above proviso is equi- 
valent to the condition p > 1. But, by hypothesis, Y > 1. Hence p is restricted 
to values satisfying 


2/2 \3/2 
[103] L<p< Qe (A-1)"-4. 


Case2. X < -—1 
If —X > Y, equation [96] becomes 


exp{ —$/2.1(—X)*°?4+2¥(—2x)'"} 


A ) 7 . r\3/2 , r\) s2 
[104] = ae 2 mpl—$v2.i( AX) +2V{2(A—X)}""): 
therefore 
[105] $4/2.{ (A—X)9/2— (—X)*?} = (2r+-))e 


where, since X < —1, 7 is any integer greater than the greatest integer (70) 
for which 
[106] $/2.{(A+1)*?2—1} > (2r+3)x. 


Let X, (r > 19) be the values of X defined by [106]. Then, approximately, 


6 mS 3/2 
lo pee \ 
[107] Y, " a Us i Sea 

2[{2(A—X,)}'— (—2X,)'"] 
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Hence when 7 is made so large that —X, > A, Y, becomes large. Remembering 
the restriction imposed upon ‘‘p”’ under Case 1 above, this shows that only a 
finite number of the terms contributed by the roots of [96] to the residue series 
for the field are significant. 
Case 3. |X+iY| <1 
We have* 
eins 1te™°0.9181(X —iV)+3(X —iY)* 


[108] FO) = eo” T+e*70.9181(X—iV)+h(X—iY) 





Hence, approximately, 
[109]  exp[jir+1.59(iX + Y) —0.73{2X Y+i(X°— Y’)}] 
i an win exp —$V2.i(A—X)°?-+2¥2(A—X)}""] 
therefore 
[110] 4V/2.(A—X)??-+1.59X —0.73(X?— VY?) = 2(K—4)z, 
K being a suitably chosen integer (vide infra) and 


tog) 8Y2A=2)""t yt 4—x)}" 


Oi ax: ea 
[111] 1.59Y—1.46X Y v1(Am) | 


This gives 





§84/2.(A—X)*” 
log (Am) 
[112] Y = Saree 
2{2(A—X)}!—1.59-+1.46X * 

The term 0.73(X?— Y?) in equation [110] is small compared with the other 
terms in this equation since P, is large and |X+7¥Y| < 1. 

We proceed tentatively, first, by neglecting this term and choosing the 
integer K so as to obtain from [110] all the values of X which are numerically 
less than unity. Taking one of these values of X and substituting it into the 
right of [112] then gives us a first approximation to Y. This value of Y and the 
value of X which has been used to determine it will now determine the term | 
0.73 (X?— Y?), to a first approximation. Including now this first approximation 
to this term in [110], this equation can now be solved again to give a better 
value of X. This better value of X now, from [112], determines a better value 
of Y, and so on. Actually the term 0.73(X?— Y?) is usually so small that it is 
accurate enough to use the first approximation to X and Y. 


3.3. Further Consideration of F; 
3.31. By [87] 


{F(pr)— F(0)} { F(pr) — F(0)} ¥2(o1 (er) ¥2(pr) { ¥2(0) 


(1 — 4i(1 — mo) | ¥2(0)}2(0k /25)o} “(= me) 


[113] F,= 


where k is as given by [81]. 
*Cf. Reference (10). 











a 
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Let 
(114] Al = ¥°{P—(Zoyo**—Z)} 
where P is as already defined* and y = 1/(1—m). The functions F are ee by 
F(p) = et [HY (3/2. (A' —X +i) *7} )/ [Ay (3V/2.(A'-X+iY)*7}] 
eins Le" "0.9181(X —iY— A’) +3 (X =i — A’) 





a8} uC Y+-e"*9.9181(% —i¥— a )+(X—i1¥—Ay 


when |A’—X+7Y| <1. When |A’—X+7Y{| > 1, then 


[116] F(p) ~ ifexp{ —$V/2.(A’ -—X+7Y)57} +7] 
if 3m < arg(A’—X+1Y) <7, 

[117] ~ texp{ —jV2.(A’ -—X+7Y)*? 7} 

if 0 < arg(A’—X+iY) < ja; 


and y2(p) is given by 

Yoo) = (fa) "2 x6 yo eft/6(4’— Xi)? . HO (R/2.(M—X+4i¥) 9} 
[118] NQ-M4 1/6 91/84’ X41 V)-" exp — hin + 3/2.(M) Xi} 
when |A’'—-X+7Y| > 1, 

95/8 
7 37re) - 
when |A’—-X+7Y| <1. 


[119] x89)! 11 —0.918le "9 A’—X+iY)} 


3.32. Approximate Evaluation of {1 —}1(1—mpo) {W2(0) }?(dk/ds) 0] 


Case (1). X K —1 
[120] y2(O) = 2-4 x6 y91/6(—- KH 47V)—4 exp{ 94/2.(—X+iV)*/*i—-hiz}, 


Ok (Am) iyo" 4 oer S/2- 
9 —) = -——|—— iV < —) 
[121] (>) “h(a X+iry*?!- 3V2.(A—X+7iY) 7}. 
But 
[122] F(O) = Na sexp{—3V/2.(A—X+i¥)*"i}; 
- 8./2.(A—X+iY)*” aa F 
therefore 
dk 9.2/9 -—-1/3, —1/3 cari /2 
7 * —2V 2.0 yo “(A—X4+1Y) “. FO) 
= —2/2.ix0 yg (A—-X+iY)'Vexp{ —3-V2.(—X+iY)7 3}; 
therefore 


a A-X+iy)"” 


*For pK pi, A! = yo3 P. 
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Case (2). X = O(1) 
¥2(0) = (39) ee yg! 9 (-X +i) er Ay /2.(-—X+iY)*"} 
_ 2s 
~ 13°°T(§) 


when |X —7Y| is small. As before 


[124] x oy! (140.9181e 79 (X—iY)} 


a = —2/ 2.0 yg 8 (A—-X+iY)'?F(0) 
ee eu - ins HY (3/2.(—X +iY)*”} 
9 aan os ua, 1/3 1/2, in /3 Se eee iae 
[125] 24/2.x (A—X+iY) ; 01279. (<X +iV)7} 
Hence 
of a 2 sei cae le 
[126] v2c0)}*( #) = —24¥2 (A—X+iY)'?.e"(-X+iY) 
. 0 « 


xX HYP 3V2.(-X+i¥)*") . HYP V/2.(-X+i¥)*} 


which is obviously large compared with unity, being of order P;'/*. Hence, in 
this case, 


\F(or) — F( OF ea) 8 | voler) v2 (or) | 


rs — $1(1—mpo) (0k/ 9s) o —mo) 
* 8 /3_.1/3 ~ 
i at monks ae 2 1 F (er) — F(0)| | Fler) — FO) 


Case (8). 1X CAH X ae 
[128] po(0) ~ 2-1/4 6 1/6 —- KX 47V)-"4 exp [34/2.(—X+7V)921—-2iz}, 
(Ok/ds)9 = —2S2.x-18 yo B(A—X47Y pe F(0) 
we — 2/2108 yg“ 8 (A—X42Y) exp{ —3V2.(—X+4+7Y)*2} 
X [1+2 exp{}-V2.(—X+71Y 2) 
ae 25/2078 yg (A — X40 ¥)* (1 —t expt —jr/2.t —8F)*}]] 





[129] = 2/2038 yo“ V8(A-—X+iY)'? 
Hence 
of ak A-X+i} | 
wont 2) Ye e exp{}v/2.(X —i¥)*” | 
X [l—iexp{ —3V/2.(X-i¥)*"}] 
[130] = -2i( aM ‘exp{f/2.(X—i¥)*"}. 


Hence, when X > 1, F; is small, being of order P,~'/? exp{ —$/2.X*”"}. 


Cases (4) and (5). X-A = O(1); X-A> 1 


Since, in these cases, 


F (pr) — F(0) 
~ilexp{ —j-V2.(Ap—X+72Y)9/2} —exp{ —j-V2.(—X+7Y)*/72}] 
[131] ~ sexp{ —jvV/2.(X —Az)*7} 
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with a similar formula for the corresponding function in pg, we have, using the 
asymptotic approximations to ~2(pr) f2(pr), 


F, = Ofexp[—3-V2.(X —Ar’)®2+ (X —Ag’)*?-+2X8]} 
and is therefore very small. 


3.4. The Factor 0; 
Using well-known asymptotic formulae for P,(—cos @) and noting that 
v; ~™ x, we find that this factor is approximately equal to 


oy \1/2 ; is : i 
[132] ( as) exp{ — i(xx OX y0 7") 0— x"? V0 7/0}. 


3.5. The Nature of the Eigenvalues 


The preceding analysis has shown that the modes of minimum attenuation 
(—Y a minimum) are characterized by the property that |¥+7Y| < 1.* It 
has appeared that the modes for which X > 1.5 (say) are very heavily damped 
by the amplitude factor F, as this becomes exponentially small like 
exp(—4/2.X*/?); while for those for which —X is large and positive, F, 
does not become small as —X is increased. The former modes correspond 
physically to rays reflected from the layer at such extreme obliquity that they 
never reach the earthf—their lack of effect upon the field series is therefore 
understandable. The latter, however, which correspond to rays reflected at 
obliquities less than limiting,f become eventually negligible at any given extra- 
optical distance through Y becoming large when —X is sufficiently increased. 

At sufficient distance, the modes of minimum attenuation dominate the 
field; since, for these, X = O(1), these modes correspond physically? to rays 
whose obliquities are near to limiting.{ Since P; is large (especially on centi- 
metric wavelengths) it appears, from the preceding analysis, that the modes 
close to the mode of minimum attenuation have a curvilinear attenuation 
factor of exp(—x!/*Y yo-?/3 8), where, approximately, 

; § 8/2. a*7\ 
[133] | i 3 logy ( 
Lyi (Am)i 


2(2A) 


i.e., an attenuation factor of 


2/3 1/3 / 3/2 
vo *.0x 84/2. A )] 
134 xp| — | ( : : 
184] exo} 2(2A)'” °8 yi(Am); 
The contribution of one of these modes to the field at long distance is, by 
[127], 
mB oi yo x! f F(pr) — F(0)} | F(pe) — FO)} 


we (a-X4+iVr)2 F(0) v2(pr) ¥2(pr) 


2x ue ‘ a y —2 fi 
x ( ) exp{ —ix—x'"(iX+Y) yo a 
wO 
i.€., approximately 
*With X as near +1 as possible. (A is large.) 


tCf. Ref. (10) (Appendices). 
{The obliquity of a ray which left the earth tangentially. 
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ae (2) vo | F(or) — F(0)}{ Fon) — FO) 
oR “* 18 


(135) (A-X+iV)” F(0) 


X 2 (pr) Yr (pr), 


where Y has the above value and X is given by [110], etc. 

For purposes of comparison the dominant term of the effective radius field 
will be given; this may be worked out separately or, alternatively, deduced 
from [98] by putting k = O therein. 

Choosing the latter method we find it to be 


[136] 7 (1—mo) valor) ¥a(pn){¥2(0)|*»Ps4(—cos 8) sec vx 


where v is the zero of ¢,_,(x)* having the smallest value of — J (v) and the 
y functions are calculated on the basis of this zero. 
Now » is given by 


¥i(0) = 0, 


i.e. by 


exp{ —{V2.(—X+7Y)*7}+7 =0 (F(O) = 0); 


therefore 


2(0) = 2M Sy (—X iV) ee et 


= = 2g th ( os gent 


8/2 
Hence the dominant term of the effective radius field is 
= eni’® ai 9 yore é 
[137] opis (2m) F a (22) (er) ¥1(eR) 


On 


4 
X exp{ —ix—x'/*(0.9272i+ 1.606) yo? 0} 


in which the values X = 0.9272, Y = 1.606 are to be used in evaluating the 
y functions. 

It will be noted that the “amplitude factor’ of our mode of minimum 
attenuation is small} compared with the corresponding factor for the effective 
radius primary mode. This is especially so when the transmitter and receiver 
are well elevated, as the rate of height gain in our case is appreciably less than 
that in the effective radius case owing to the smallness of our value of Y above 
in comparison with the value of Y (1.606) appropriate to the first mode of the 
effective radius field. This, however, is more than made up for at long distance 
by the comparatively slow rate at which our field decreases as the range 
increases. 

It is also interesting to note, from the above, that the field at long distance 
varies exceedingly slowly with wavelength (compared with the rate of vari- 
ation of the effective radius field with wavelength) and shows a tendency to 


‘ 


*Here, the ‘‘R" in x = 2xR/) is to be the effective radius. 
{Pe T. 
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increase slowly as the wavelength is increased. This is in accordance with 
what is observed. The field also varies only slowly with (Am). 

At relatively short ranges (i.e., at those ranges where a field of the order of 
magnitude of the effective radius field becomes important), the individual 
contribution of these modes of minimum attenuation* is comparatively small, 
owing to the smallness of their amplitude factor (vide supra). However, at 
such short ranges, great numbers of the eigenvalues for which X = O(1) and 
for which X < —1 become important, owing to the sluggish manner in which 
Y increases with |X|. Summation of these modes at such short ranges would be 


TABLE II 
Rate of attenuation of ob 
received power (field 
strength squared) at long 
distance 
Inversion (decibels /mile) 
thickness —_ + - 

h (ft.) (meters) (Am), \ = 6 meters » = 10 cm. 
7,623 100 2 0.181 0.308 
25 8 0.136 0.264 
10 20 0.106 0.234 
11,4385 100 2 0.16 0.281 
-25 8 0.13 0.242 
10 20 0.10 0.216 
15,246 100 2 0.15 0.251 
25 8 0.12 0.218 
10 20 0.10 0.196 
5 40 0.180 
19,058 100 2 0.15 0.234 
25 8 0.12 0.204 
10 20 0.098 0.186 
22,870 100 2 0.145 0.228 
25 8 0.113 0.200 
10 20 0.094 0.180 
5 40 0.166 
26,682 100 2 0.136 0.215 
25 8 0.110 0.188 
10 20 0.094 0.171 
5 40 0.158 
30,494 100 2 0.12 0.200 
25 8 0.10 0.176 
10 20 0.08 0.156 
5 40 0.146 





Note: (1) The values obtained by Carroll for his model are 0.15 db./km. (0.24 db./mile) for 
= 10 cm. and 0.10 db./km. (0.16 db./mile) for ¥ = 6 meters. 

(2) The last height in the above table is the height of the du/dh discontinuity occurring in 
Carroll's model. If we use (Am) = } with this in the preceding theory we have a layer whose param- 
eters are equivalent to those at Carroll's du/dh discontinuity. Calculation of this case yields for 
the far-field received power — distance attenuation the values 0.24 db./mile for = 10cm. and 0.17 
db./mile for \ = 6 meters, agreeing almost exactly with Carroll's figures. This illustrates the point 
previously made—that his model is really equivalent to a sharp elevated layer. 


*We mean, of course, the modes near enough to that of minimum attenuation to be important 
and the mode of minimum attenuation itself. 
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very tedious and as Carroll has already carried out a similar kind c. calculation 
(3), verifying that this mode sum does indeed give an exceedingly close approxi- 
mation to the “effective radius” field at short range, there seems little point 
in carrying through this calculation here. It is almost certain that we should 
likewise find a close fit with ‘‘effective radius” theory at short range as our 
equation for the eigenvalues is of exactly the same type as that encountered by 
Carroll. 

The rate of curvilinear attenuation [133] decreases when the layer height is 
increased and also when (Am) is increased. In the cases likely to be of impor- 
tance in practice the layer height is a good deal less than the height of the 
u’-discontinuity in Carroll’s model but (Am), will be appreciably greater. 


3.6. Numerical Results 
A few values for the rate at which the mode of minimum attenuation is 
attenuated have been worked out and are exhibited in Table II. 


4. CONNECTION WITH RAY THEORY 


4.1. Let 6, be the greatest angular separation of a surface transmitter T from 
a surface receiver X consistent with a singly reflected ray from T to R. Following 
the terminology introduced in the writer’s previous paper (10), we shall refer 
to the distance over the earth’s surface defined by 6, as the “limiting distance’. 

If the atmosphere below the layer were a vacuum then the value of 6, 
would be 2(24/R)'/?. However, we are assuming that the refractivity profile 
below the inversion layer closely approximates to the geometrical theory appro- 
priate to a vacuum above an “‘effective’’ earth of radius R’ = Ryo. The value 
of 6, to be used is now 2(2h/R’)'/? which gives a limiting distance of 2(2/R’)!/2, 
i.e., \/Yo times the limiting distance over an airless earth. 

A ray reflected once from the inversion will in general return to the earth 
over an angular distance 0, where 6 < 6,. Let ¢ be the actual angle which the 
ray makes with the inversion layer; by simple geometry 

_ ©08($— 36) 


h 


cos @ 
i.e., since all angles concerned are small, 
[139] 16,2 = 6¢6—16° 
giving 
[140] 6 = 2{o— (?—10,2)/2}. 


The general wave equation satisfied by the field functions is, in two dimen- 
sions (Oy horizontal, Oz vertical), 


Cit ee 
[141] ay’ aett 2 i = 0. 


If x varies sufficiently slowly with height, 
[142] x = xoexp{ (tp/C)x/?(y cos ¢’+2 sin ¢’)} 


is a closely approximate solution. This represents a local ray (whose direction 
is that of the normal to this plane wave) which is making an angle ¢’ with the 
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local horizontal direction.* Expressing x in terms of a height function f(z) 
(to parallel our process of separating the variables r, 6 in our wave treatment), 
we write 


[143] x = f(z) . exp{ip/C) x”? y cos ¢’} 


so, to this approximation (neglecting the derivatives of x), 






2 
(2) + a x sin’g’. f(z) = 0, 





i.e. 


[144] f"(2) + a ¢” f(z) =0 





since, in the propagation problem, ¢’ is small and x is near to unity. Comparing 
this equation with [18] we see that, to a particular mode (particular permitted 
value of «), there corresponds a ray whose local angle of elevation ¢’ is given by 
[145] g'? = 2(p—e,) +6 


where e, = #(e) and .Y (e) is omitted because it is always small compared with 
2(p—e,)+6 for all modes which are significant (§7.32, Case 3).f 

This gives 
[146] do" 








—2(€,— 360) 


= —2Qyx-2/3 y,2/8 










oe 







X, denoting the value of X (see [96]) for a particular mode order r. 

Thus only modes for which X, is negative can represent rays which could 
have originated from a surface transmitter. However, for an elevated trans- 
mitter at height hy (p = pr, P = Pr), the value of ¢’ at the transmitter, ¢7, 


say, is given by 








gr? = 2(pr—e) +67 
= Qx-2/3(P,— 92/3 X,) 






[147] 
since, for all moderate transmitter heights, hy K Réo, giving 67 = do. Let 
[148] Py’ = yo°*8 Pr, 






then only those modes for which 
[149] ae SPS 





have associated rays which could have originated from the transmitter. Not 
all of these rays will, however, necessarily reach the receiver; for them to do 
so we require also (with an obvious notation) 

[149A ] Ae SP. 







Remembering that, for all significant modes, X, is certainly less than 1.5 
we see that the condition for all rays associated with significant modes to be 
capable of travelling from T to R is that 
[150] min(P,’, Pp’) > 1.5. 


*We reserve @ for the value of $' at the inversion. ; 
+[145], incidentally, agrees with the law of variation of ¢' with height required by Snell's law, 


for pr sin $' = poo Sin do is equivalent to g'? = go? +2h/R+ (5 —So). 
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If, however, this condition is not fulfilled, the rays corresponding to values of 


X, given by 
{151] LS Xx: = max(P7’, Pr’) 


do not reach the receiver but pass near it.* These modes must therefore con- 
tribute to the received field through diffraction. However, there are always 
important modes (i.e. given by |X| < 1) whose rays do reach the receiver 
direct}, for instance those for which 

[152] —1 < X, < min(P7’, Pp’). 


Thus in this theory, height gain is to be explained by the fact that it is 
possible for the rays of every significant mode to travel from transmitter to 
receiver only when the heights of both make [150] true, i.e., when the heights 
of both exceed 
[153] 0.44 R’¥* 24, 


4.2. On the Relative Phases of the Modes 
We have 
6, = 2(2h/R’)'”? 


[154] yo¥? 26/2, PU? x3, 


so that 
(155] P¥24/2 = Ayo!/2 x18 Oy. 


Since P and A are always large, the separation of successive eigenvalues 
(X —iY) is small and, when the ‘‘angle of fire’ @ is not too small (i.e., when 
—X,> 1), we find from [105] that the difference of the real parts X of adjacent 


modes is 
T 1 
(156 V3 KR 


Tr 


In terms of ¢ and 6, this is equal to (cf. [95] and [145]) 
T 1 
ii ay = (—) 
Hence the actual phase differencet between adjacent modes is 
yo cp eseieeieerees 
¢,—(¢, — 701 ) ‘o 
[158] = 2176/0, 


by [140], where 6, is the angular separation between points on the earth’s 
surface struck by a ray for which the mode number is r (r is as defined in equa- 
tion [105]) and @ is in this formula the angular separation of transmitter and 


*Because max (|Pr’ —X,|; |Pr’' —X,-|) = O(1). 

tI.e. after one or multiple reflections. 

{The phase difference between the amplitude factors F, is of the order of 5X, itself and so can 
be neglected in comparison with [158] when x" #9 is large, i.e. at distances well beyond optical. 
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receiver measured relative to an “effective” earth, i.e. R’@ = (transmission 


distance). 
We show immediately that this formula holds also for small ‘‘angles of fire”’ 


oo (when |X,| < 1). By [110] 


6X = JB vice 
~ 9/2 (A—X,)'”—0.56 
T 1 


V2 (A-X)" = (-X) 


if \/A is not too small,* i.e., when ray theory is applicable, 


cs 1 
[159] ~ 77,17 6 
as before. 
The form of [158] shows that we may, in computing the phase separation of 
successive modes, measure both @ and 6, relative to an actual earth. 
Note.—In the above, we have assumed that the transmitter and receiver 
are both on the surface. Since hz and hg are both certainly very small compared 
with inversion height h, the above results are applicable to a close approxima- 
tion. However, since it is possible for rays proceeding from an elevated trans- 
mitter J to an elevated receiver R to never reach the earth’s surface, it is 
clear that, for such rays, 6, has to be defined with reference to a sphere through 
R concentric with the earth. In such cases 6, is very nearly 6, and, in fact, by 
carrying through the above analysis again for these kind of rays it is easy to 
see that the phase difference will work out to 276/8,. 
4.3. The Difference between the Attenuation of Successive Modes 
By [107], the attenuation in nepers over an angular distance @ when 
—X,> 1, i.e., when x18 = 2 19,7) > 1,f is 
f 4x¢,° \ 
00] 510g Sl (Am), } 


where @ is the angular range measured relative to an “effective” earth and is, 
of course, the angle at which the ray corresponding to X, strikes the layer. 

But when |X,| < 1 we have, approximately (since A is large for the inversion 
heights and wavelengths in which we are interested), 


Axor 
Te X,))72 8 log S28} 








J lg| *} 
~ F178, 78%, 8 Cam) S 
[161] > aie (tee 


since in this case 6, = 0,, 2¢, = @,. Hence the approximate formula [160] is 
generally applicable. 


*Since —X, <1 here. 
$8 yal 8(p? 20,2919 m 3/8 49 9(G, — 467) = 219 70"? rceragay > 1 dy [146]. 








334 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


The change in the attenuation [160] when the mode order is changed by 
unity 

om 4xg,*\ 
162 = = logy on 
[162] log 
since the logarithm is usually moderately large and ¢, changes relatively 
slightly when r is increased by unity. 

Now, when X, < —1, we have 

xy 9\? 60, = 8(2(A—X,)'?—(—2X,)'"} 





ae oe oe 
“hua” Gea 
i ee 
~ 08x og, (67-2) 
: ee ree ae 1 
[163] vi e "ee" (or —7Or° —9o,")'” , 


Therefore the actual difference in attenuation between the adjacent modes 


- 6 on} Axo, \ 


(Am), 


xY0 or (o; —16, ye 


ge 8 <br 9 
[164] - ogee {eee 


This is, of course, only true when —X, > 1,* i.e., when @, is not too near @;, so 
this attenuation is actually less than about 








s 2a fot 
[165] xy 8 6,> 9 ton (Am) 


Otherwise (|X,| < 1) we have, by [112], 
log) 8¥2-(A=X)*"\ 


3 = (Am) f 1 \ - 
os TIC n ee oats ox 


84/2.(! or 
log SX (Am) 0.730 . 
Or x eye —— xBy'8 6 b,— (67-107) | 


Hence in this case the actual difference in attenuation between the modes is 


(1460 a geen 
vs yo* $,°0 (Am) 


_ 5.847 86 xOy ce 
- = rn el 


Hence this attenuation difference is at any rate less than about (2/P) (6/6) 





ll 





1/3 yo!/3(g,2 — 40,2) 2m x3 x0 V3 A(% -0,) Sa: 
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times the logarithm. Since P is so large this shows that, for all practicable 
ranges, the magnitudes of many modes near to any given mode must be 
practically the same. 


4.4. The General Characteristics of the Field at Long Distance 

4.41. It will be convenient to interpolate at this point a word about the 
classification of the various field regions. When the layer is absent the field at 
points out of optical ‘‘sight’’* of the transmitter is due solely to diffraction from 
rays outside, and upon, the tangent cone drawn to the earth’s surface from the 
transmitter. Such points may properly be said to constitute the ‘‘diffraction” 
region of the transmitter. When, however, an elevated layer is present energy 
can enter this region via rays reflected once (or more times) from the inversion 
layer. We shall therefore speak of this region as a whole as the ‘‘semidiffraction”’ 
region and further subdivide it as follows: the ‘‘primary”’ semidiffraction region 
is that part of the semidiffraction region which can be reached by rays reflected 
once only; the “‘secondary”’ semidiffraction region is that part which can be 
reached by rays reflected twice only, and so on. 

Note that the distance range.up to the end of the primary semidiffraction 
region is what we have already defined as the “‘limiting’’ distance. It would 
now be more logical to call this the “‘first order limiting distance”’; the distance 
up to the end of the secondary semidiffraction region as the “second order 
limiting distance’ and so on. 

The factor depending on the angular distance @ which is contained in the 
expression for the field contribution of a given mode is 


vim | 8 {sx6,"\ ; | 
6 exp] 6, 8°. (am) +i, |, 








1.e. 

2 ~—1/2, ) a i 6 ys Oat A,” y'} 
[167] 6 exp) Vr 6, lal 57 (2 +6, f 
where the phase y, equals x!/3 yo~?/3 0X, 

= 1-1; 6, , 

[167A ]f = ~39 XxYo o( 9-0.) 


and, by [166], the phase separation ¥,41—y, between adjacent modes is 276/8,. 
This formula is validf for 6, > 36,°, i.e. for d,; > 24h where d, is the total 
distance covered by a single hop of the ray ‘‘0,’’. 


4.42. The Field Variation within the Primary Semidiffraction Region 

As this region is entered we would expect the field strength to drop relatively 
rapidly, owing to the loss of direct transmitter-to-receiver rays. Furthermore, 
since the effect of rays which have been reflected from the elevated layer is 
extremely small at short ranges of penetration into this region (owing to the 
largeness of the angle of attack of such rays upon the layer) we might reason- 

*Over an “‘effective”’ earth. 


tA ppendix 9. 
tA ppendix 10. 
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ably expect that the law of field variation at such short penetration will follow 
that of simple diffraction around an earth of ‘‘equivalent’’ radius R’. As has 
already been remarked this has been verified to be the case by Carroll for the 
bilinear model* of refractivity profile considered by him. It is interesting to 
note that the present theory also indicates that there will be a sharp drop in 
field as we enter this region although it would, of course, need detailed and 
tedious calculation of the sum of many modes to verify that this drop did, in 
fact, fit the ‘‘effective radius’ theory curve at these ranges. For, if 0/6, is 
small enough, practically all the modes significant in the field strength calcula- 
tion are (to all intents and purposes) in phase, as their phase separation number 
2176/8, can be made as small as we please by sufficiently decreasing 6. When, 
however, the range is increased so that the receiver moves beyond optical range 
(or, in the case of the transmitter T near the surface, into its diffraction field 
of classical theory), @ becomes certainly as large as O(x~?/*), making 276/06, 
larger than O(P-"/?), since 6, < 6,. For ranges which represent a relatively 
small penetration into the primary semidiffraction region, the number of modes 
giving significant contribution to the field is large, owing to the relative small- 
ness of 6. Careful consideration of theory already given shows that this number 
is of order P. Hence the total phase shift over all these modes equals O(P'””), 
at such distances, and this is quite large compared with 27. Hence the very 
large group of modes which were effectively ‘‘in phase’’ within optical range 
will, shortly beyond this, contain a large percentage of interfering elements. 
We thus have a qualitative explanation of the phenomenon in terms of the 
present theory. 

As 6 is now increased beyond the relatively small penetration of the partly 
lit regiont given by @ < 36,? the effect of the rays reflected from the inversion 
begins to be important. We note that rays for which 6, = @ (rays proceeding 
from T to R via a single reflection) are exactly in phase, as are also the double- 
hop rays (6, = 36), the triple-hop rays (6, = 36), and so on. Owing, however, 
to the high angle of attack of the multi-hop rays their contribution to the 
received field is unimportant compared with the contribution of the single-hop 
ray. For example, the ratio of the field strength due to the double-hop ray 
alone to the field due to the single-hop ray alone is (by [167]) approximately 

1 (Am) 6° 


[168] 2/2 “-, 


which is <1. 

As we have seen, the total numberf of admitted rays (modes) is very large, 
especially for the shorter wavelengths. Accordingly, the dominant part of the 
field due to the various modes of spans 6,, at a given angular range 8, can be 
expected to come from the mode group or groups which make 276/68, as near 
zero or a multiple of 27 as possible. Thus, at points not too near the beginning 


*As has already been remarked his model behaves like an elevated layer. 

{The primary semidiffraction region. 

{The number of terms in a mode sum adequate for a closely approximate representation of the 
field in the region under discussion is at least O(P). Further out, of course, this number decreases, 
the are modes clustering ever closer around the mode of minimum attenuation as range is 
increased. 
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of the primary semidiffraction region, the modes represented by rays close to 
the geometrical “‘single-hop”’ path are dominant. 

At points near the beginning of this region {@ = O(36,?)} the attenuation 
of modes near that of minimum attenuation is approximately 


8 acts!) 

(A) 8, lel Am _) Mepers 

while that of modes near the single-hop mode is approximately 
xp ) 

(B) to( 52 nepers. 


Both groups of modes are effectively in phase (since 2170/6, = O(@,) K 1 
for modes near that of minimum attenuation). Hence, since expression (B) 
is obviously appreciably greater than the expression (A) it follows that, for 
small penetrations of the primary semidiffraction region, the many modes near 
the mode of minimum attenuation dominate the contribution from modes 
near the single-hop mode. Hence, except near the beginning oi the primary 
semidiffraction region, the field is dominated by the group of modes whose 
representative rays are near the single-hop path from T to R. Since for such a 
mode group, 6, = @, its law of variation with range would appear to be described 
by the function 


= f x (% ’ ” 
[169] gx \Stam) “+0 \ 


It therefore appears that the field increases slowly with 6 for 6 less than a 
certain submultiple of 6, (but, of course, @ must be >36,? for this law to hold) 
and decreases when @ exceeds this multiple of 6; (@ < @,). Actually, we have 
omitted the ‘‘amplitude’’ factor from [169]; this is allowed for in Appendix 11 
and does not affect the phenomenon we have just remarked upon. We thus 
conclude that, when @ increases in the range 
[170] 30,7 < 6 < 0.58726,, 


the field strength - range curve first flattens out from the initial steep ‘‘effective 
radius’’ descent and then begins to rise slowly until 6 = 0.5872 6,. Thereafter 
the field gently decreases to the end of the primary semidiffraction region 
(@ = 6,). This can be explained physically as follows: when we first reach 
the region where the modes near the single-hop mode are dominant the 
incipient field strength increase is explained by the fact that the angle of 
attack ¢ of the single-hop ray upon the layer is decreasing so rapidly that the 
consequent rate of increase of its reflection coefficient more than offsets the 
steadily decreasing ‘‘free space’’ term. Near 6 = 6,, however, d¢/d6 approaches 
zero and the field then slowly decreases. 


4.43. The Field beyond the Primary Semidiffraction Region 
The secondary semidiffraction region (0, < 6 < 26,) 


At the commencement of this region the most powerful modes (6, = @;) 
are in phase (so also are the multi-hop modes 6, = 36z, 30z,... but, as we 
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have seen, contribution from these is negligible compared with that from 
6, — 6,7). 
As @ is increased past 6;, the attenuation of the modes of least attenuation 
(for which 0, = @,), viz. 
0, JS 4x0,’ + 
171 log 
en Or s\ (Am) 


increases steadily, and so also does the phase separation between adjacent 
modes in this group, 27(6/0,—1). Thus the actual rate of loss of field strength 
in the first part of the secondary semidiffraction region may be expected to be 
a good deal more rapid than the rate of loss of the modes of minimum attenu- 
ation alone. This relatively rapid decay in field strength may be expected to 
continue until the phase separation of the adjacent modes begins to decrease 
with 6, i.e., out to @ = 36,. At this distance the modes which are in phase are 
and so the attenuation of the field at this distance is roughly 
~ 3 

[172] 2 log) * Suit log (3 01). 

As @ is increased past 36, and on towards 26,, the modes of minimum attenu- 
ation (8, = 6,) become more and more nearly in phase, reaching complete 
phase agreement at @ = 20,. At @ = 26,, therefore, the attenuation is approxi- 
mately 


given by 6, = 561 





‘ 4x0," 
8. (am) ft? 


i.e. very slightly less than the attenuation 6 = 30;. 


[173] 2 log + log (26z), 


The more distant semidiffraction regions 
Following similar reasoning, the field decreases rapidly in 20; < 6 < $6,, 
remains practically constant in $6, < 6 < 36;, and so on. 


4.5. Note on the Fading of the Signal 
1.51. On this theory, fluctuations in time of the received signal will be caused 

by moment-to-moment variations in the inversion parameters (Am) and h. 
The relative variation in each of these quantities can be expected to be small, 
for an inversion which is in a reasonably stable condition. Since (Am) enters* 
only through the imaginary part Y of the eigenvalue the effect on the distant 
field of small variations in it is relatively slight. This, however, may not always 
be the case as far as variations in h are concerned because the relative phases 
of the modes are also affected. We now proceed to study the effect of small 
variations in h. 

(1) Within the primary semidiffraction region 

Except near the beginning of this region it is the modes near the single-hop 
ray 6, = @ which are predominant. When hk is changed these modes still exist 
and are still in phase so the change in field strength at @ is, according to this 
theory, due only to the change produced in the angle of attack of the single-hop 


*This is not strictly true but it is a very good approximation. 
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ray upon the layer by the change in layer height. For 6h/h small, this change in 
field strength is extremely slight (as may be verified from [169], remembering 
that 4(éh/h) = 60,/0,). 

Hence within the primary semidiffraction region there should be very little 
fading. 

(2) Outside the primary semidiffraction region 

Here, the modes near the mode of minimum attenuation are the important 
ones. The increase 6h in layer height produces a number of extra modes of this 
type and, by Appendix 8, this number is approximately 0.9P!/? .6P when P 
is large. The increase in relative phase separation (applicable to these new 


modes) is 
- ast) _ _276 
[174] i( 6). ~ 6: 661, 
9 
= ae 60, atdO= 02, 
L 


a 66, atéd= 5 Op. 
Or 


For noticeable fading we require that the total phase spread of these new 
modes should pass through at least one multiple of 7 as the laver height vari- 
ation 6/ increases from zero to a magnitude physically plausible during the 
period of the fade. Reserving now the symbol 6h for this latter magnitude we 
have, as the condition for noticeable fading, 


te = xX P'" 6P > say, 
Or 
Le. 
[175 6h > 0.4 R42 h/4 = (say). 


A simple calculation now shows that condition [175] can quite likely be 
satisfied by 5h for centimetric waves but that this is not at all likely for meter 
wavelengths. 

Hence we would expect appreciable fading beyond the primary semidiffrac- 
tion region on centimetric waves but hardly any on meter waves. 


Note.—In the above we assume that the /ayer itself is in a reasonably stable 
structural state. Also, we neglect the disturbing (random) effects produced by 
atmospheric inhomogeneities in the path of the downcoming (reflected) 
wave—these are bound to produce time variations of signal strength also. The 
above theory must therefore be applied with caution. 


4.52. Limitations of the ‘‘Layer’’ Theory 

If the layer thickness becomes so small that it is comparable with the local 
scale of turbulence it will no longer be possible to treat the layer as horizontally 
stratified, nor, indeed, as a statical entity at all. In this case the layer could 
probably best be thought of as a region where turbulence was particularly 
intense, i.e., we should have what is called by Gordon a ‘‘turbulent layer’. 
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The effect of such a layer on radio propagation has been studied by him (7) 
and it seems likely that it could, in fact, give substantial contributions to the 
long distance weak fields. In this connection it might be helpful if the scattering 
theory could be extended to take account of contributions from the turbulent 
layer above the mid-point of the path when this mid-point is out of optical 
“sight”’ of the transmitter or receiver, or both. Such considerations are, how- 
ever, beyond the scope of the present paper. 

The preceding theory also assumes that the layer is a spherical surface con- 
centric with the earth. In practice this requirement is not, of course, satisfied, 
but (cf. Reference 10) the theory can be expected to apply if the layer is well 
developed above the mid-point of the path with horizontal extent large com- 
pared with the linear dimension 0.1 R*°/*\?/*/h'/? where h is the layer height, R 
the earth’s radius, and \ the wavelength of the radio waves. 


4.6. Note on the Possibility of Identifying the Principal Agency Responsible 
for the Distant Fields 

(1) It has been noted that, if the field is due to elevated layers, then there 
is a second region of rapid field strength decline shortly beyond the limiting 
distance 2(2hR’)'/? (the first such region is the ‘‘effective radius field region” 
which begins just beyond optical range). It is not likely that layers of the kind 
we have been discussing will occur higher than about 30,000 ft.; hence, if 
there is no such second region of rapid field decline within about 480 miles, we 
may safely say that the phenomenon is not due to an elevated subsidence 
inversion layer. If, however, such a second region occurs, the distance of its 
onset from the transmitter will enable us to obtain a fair idea of the height of 
the layer responsible for the propagation. 

(2) The observed field is Rayleigh distributed with time. There is no reason 
why the layer parameters should have to vary with time so as to produce this 
effect. 

(3) Signals on 10 cm. waves personally observed by the writer at Toronto, 
emanating from a transmitter at Ottawa, seemed to show much too deep and 
rapid fluctuation with time to be explicable in terms of time variations of 
parameters associable with any physically plausible layer. 

5. CONCLUSION 


It has been found that the ‘‘standard atmosphere’, which is characterized 
by a steady and slow decrease of (u—1) with height, cannot by itself give rise 
to the long distance V.H.F. fields observed in the deep shadow region of a 
V.H.F. transmitter. The analysis of Part II has, however, shown that high 
level inversion layers, of strengths and heights which could occur in practice, 
are capable of producing the effect in question (cf. Table II), but, in view of 
§4.6 (2), it seems almost certain that they cannot be the sole agency respon- 
sible, (2) being characteristic of the field received through scattering from 
atmospheric turbulence. It therefore seems probable that the ‘‘usual’’ signal 
levels observed must be due mainly to electromagnetic scattering from 
atmospheric turbulence, unusually strong* signals being caused by high level 
inversions. 


*I.e. signals whose mean over a period of hours is large compared with the signal mean over a 
period of months. 
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APPENDICES 


Appendix 1. JS (e) 


The equation satisfied by u,_, (é) is uy _y+x?[2(p—e) +5}u,_, = 0, the dash 
denoting differentiation with respect to p. Hence, since p is real, 


vy +x" {2(p—2)+65}05- == (): 


therefore 
U5 -yU5—4 —V5—-4Us_y = 4x"i_g (e) X Us—iV5-4, 


and 


, ’ es 2- i. 2 
Ms—sVs-y— UsMy | = Av (€) lies dp. 
p= 


Put s = », then Uy, vanishes at p = 0; so also does v5-, (at p = 0). Also 


u(t) ~ e~*, v(t) ~ ef, 
as § > + ~. Therefore 
IGY a ome gg 
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Appendix 2. Langer’s Solution 

This solution reduces asymptotically to the W.K.B. solution for large values 
of the argument; it is finite at the “turning point” at which Q = 0, and reduces 
to form [20] near the surface and [23] at large heights. 


It is 
U 7 
( ae (W). 


Q’=W 


Now 


so that, taking 


W 1/2 ; 
tin (e) J,(W) 


and using Bessel’s ksi for a function of order 4, we find 


gee * c at we en. 
Q” = —(W’) oS aant +viw)( 415 Gy? ap pAsev 
Sai ey a) “| — 1w'\, 
ae “SA Ws TAPS oh 


tw 3 2 oe 18 
-(o- 361 14 OF -12)o. 


Appendix 3 
dF_1(, db dys) 
dp rte dp Ms dp 
im = SPW)! Hy? (W)}'— HYP (W) (Hy? (W)}] ow 
2ix 


vo 


Appendix 4 
The validity of §4 depends on condition [424A ] being satisfied. This condition, 
however, is approximately equivalent to 
4/3 
Mo 


|jp>—e+36| < 1.5 (say). 
0 


Since .Y (e€) = O(x7?/%) this is, effectively, 


4/3 
Mo ‘ 
; UP 


The left-hand side has a maximum where 
m(1—m)+3m'(p— .# (e) +36) = 0. 


But 
= O(myo?/50?), 
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therefore 





4/3 


m we 3 
max 5 te # (€) +36} = O(m'"*) (5-8. sme) 
0 






<1: 










Note.—For the exponential profile the maximizing value of p is approxi- 
mately* 






3 50 
8 Mo 








so that 


4/3 
m 3 1/3 1 
max] "to— 2) +48) | = 5M € 
0 





s2 





e 1, 
<= —— $5 


8e 
Appendix 5. Condition for the Solution to Appreximate to the Watson Airless 
Earth Case 
We have (for the eepenetia weg 


os _<M pr - Qmp* Leo 
J, itm l—m)') F 
1 36 ») f 2mo \ Q? 
a Ze 1 smna?+ 25 m exp) "7, Py. F, dP. 


We have to approximate to this when Z» — 0; noting, however, that 
2(0) = 0, 
2’(0) = F’(0), 













this integral approximates to 


oe WS ,) =i 5 
120 ds Ses TS 


= $2) X Fp(0). 














= 2V2.(-E+1Z0)4 exp! 
= —2/2.(—E+}Z))*i 
— 2/2. (99 /8+/2)). et? | 
The integral approximates to 
3(9r/2)? a 2. 


Therefore the condition for the propagation to approximate to “‘airless 
earth” is 





—19/2.i(—E+42Z,)*”} 








Zo<i (say). 


(Note that this condition is independent of mo, which is what we should expect.) 





* For this p— Bie) +}b6— p. 
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Now the condition for propagation to approximate to ‘‘effective radius”’ is 
my?/Zo <4 (say), 
1.e., 
Zo>%} 
if my) = 3, as in the “standard atmosphere’”’ case. 

Hence, for a standard atmosphere, there is no value of Zp for which the 
propagation approximates neither to “effective radius” nor to ‘‘airless earth’’. 
Hence, for the particular type of profile with which we are concerned in §4, 
propagation approximates either to ‘‘airless earth’ (Zo < 4) or “effective 
radius” (Zo > 4). 

Appendix 6 

If + is the measure of thickness of the layer (thickness through which m’ is 
approximately constant) and is large, then, bearing in mind that the smallest 
value of (p1;— +48) is roughly of order x~?/*, condition [65] for the layer to act 
as ‘‘sharp”’ is, effectively, 

7K RY 23(Am), 
or (Am), > 7/R'3 n2/, 


It seems likely that this condition will often, in practice, be satisfied. 


Appendix 7 
Suppose p < 1, then 
aici Ce Bae aut ss Zz 3/2 
F(0) = 8/2. x(pi—e + 481 sda exp{ 3V 2.X1(p1 e+ $61) 1} 


_ (Am): (Am) mpl v2. iinet He) (oye oh) lt 


x add not iin)" 
p2—€+ 352 


X exp{ —3V/2.xi(o2—e+ 462) °?y2} 


m | exp{ —$/2.2i (ore 48s)” 
= go 3V2.X1 (pe e+ 362) 2} Se (Am), 


a a) aging SF 
2 2%: 


showing, since ao > 0, that F(O) is very large indeed. But, analysis 
similar to §3.2 will show, from (A), that the modes which dominate the field 
at long distance make x?/*(¢— 369) of order unity, which in turn makes F(0) of 
order unity. It is therefore inadmissible that p should be <1. 


Appendix 8. Approximate Solution of Eigenvalue Equations 
Case |X| < 1 
We make use of the fact that A is large. Then [110] is approximately 
0.3001 A*/?— (0.4502 A'/?—0.2531)X = K—0.3; 
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(0.3001 A*’?+-0.3333) —K 





therefore ~ "0.4502 A"? 0.2531’ 
i te 
me a 
1/2 

2@a)'*—1.59+| 146-( 2 ) kx 
3/2 1/2 
og 82-4} a 4s-(2) | 
= 5@4)*—150 |!~Ziogcs) UI -a@my*—150 





aoe ( 2 T.. 


“sae et Ga 1.59 0.4502 A? 0.2531 





K being restricted to make |X| < 1. 


Appendix 9 


—2/3 2 
or (surface) 


= 4x79? 8 (g,—46,)” 


— 1,2/3, -2/3 te .) 
ax yo -6\6, 9 


eae 1 wee: %, . 
x1 OX, = — 35 x0 v%9(82"_o,) 


r 2/3 
—-X, = by? yg 


by [139]. Hence 


the ‘‘0’’ here being measured (as in the text) relative to an ‘‘effective’’ earth. 


Appendix 10 
The ray theory discussed in the text rests upon the assumption that the 
angle of attack of the general r™ ray upon the layer is small enough for us to 
(i) replace the exact equation [138] by the approximate equation [139]; 
(ii) replace sin ¢’ by ¢’ in [144]. 
This requires that 
dr <7 say, 


i.e. that 





Be iets on 
49,82 +9) < 45) 


i.e. (effectively) that 


6, > 367”. 
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Appendix 11 
We wish to find the turning point of the function* of @ 


2 —3 2\ 2) —-1/2 
rtf) fur min E8) 


2\2)-—1/2 
tig yg ele eis eae +801) \ . 
L 


i.e. of 


Taking my = } the en differential is 


__ 60 6,'6(6° +01") 
ase pet? ea HO —eey THOTT 


If, as a first approximation, we neglect 6?/6,? as small, we find that the value 
of 6 giving the turning point is 0,/+/3. Thus this first approximation is justified 
and the next approximation gives the value of @ for a turning point as 0.5872 6,. 


*N.B. by [123] F contains the factor 


X,+iY,\ 2 
146 — my (ARE t 


= }1+(1—mo)*( A=) ‘~ 


=fita- seal CG) ‘oe 


and 6, = @ for the single-hop mode. 


of which the modulus 
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III. THE CASE OF TWO ELEVATED LAYERS! 


By Francis H. NORTHOVER 


ABSTRACT 


The general theory developed in Part I and in the beginning of Part II is 
applied to derive an equation for the eigenvalues of the propagation for the 
case of two sharp elevated layers. This equation is very briefly discussed and 
some special cases noted. 


The solution we have developed in Part I* (see equation [67], Part I1*), being 
expressed as it is, in terms of a height integral, can readily be extended, if 
desired, to cover cases where there are any (finite) number of layers present. 
It is, however, considered relatively unlikely that more than two layers would 
be simultaneously present, so we shall, for simplicity, confine ourselves to 
this case only. Let the heights of upper and lower layer be, respectively, Rp: 


and Rpo. 
Then, for p < p; but >p2, we have, by the analysis of Part II, §1.1, 


Q(p) = Po)+ saa exp{ —3V2.xi(o1—e+ 461) v1} 
[1] = 1(p), say. 


For the important modes it transpires that t@ {x”’*(e—}50)} is of order unity 


or <—1; we may therefore use the approximate formula 
[2] F(p) = iexp{—3V2.(o—e+48)*"ixy} 


in the above. 
Applying the integral formula [67] and using the factf that the contribution 


of the ‘‘standard atmosphere” portions between the layers, to this integral, is 
negligible in comparison with the contribution of the layers themselves, we 


have the approximate equation 


" — Sia” exp{3V2.xi(p2—e+ 452)” 2] {Qy (p2)}" 
+ exp{ —4/2.xi(p1—e+ 481)? ?91}. 
Let 
3/2 
re ihe 2.xi{ (pr—e+ $1)? 71— (p2— e+ 452) *¥2} J, 


so that, when p > 1, 
1Manuscript received February 1, 19565. j ; : 
Contribution from the Department of Mathematics, Memorial University, St. John's, New- 


foundland. 
*Part I.—Can. J. Phys. 83: 241. 1955. Part I11.—Can. J. Phys. 33: 316. 1958. 


tPart II. 
{Part I. 


347 





348 CANADIAN JOURNAL OF PHYSICS. VOL. 33 
° ‘tes 3/2 
[4] Qi(p2) = 7 exp[3v/2.ix(p2—e+ 362) "/"72] 


and, when p <1, 


. (Am) 4 ‘ 3/2 
(1) = s355p s exp[—$/2.xi(p1—e +4 
[5] bi (1) 8./2.x(pi—e+ 48)” exp[—3-V/2.x1(p1—e+ 361) Yi] 

It is only possible to reach clear-cut resuits when p > 1, or when p < 1. 
It is shown in Appendix 7 (following Part I1) that the case p < 1 cannot occur. 
We therefore consider the case p > 1. 

In this case the equation for the modes is 


; (Am), eed 3/2 
os Md i. <4 yiz expl— aV/2.xi(pi—et+ $61) 14 
(Am)> 3/ 
ie ctaisins ee a 15, /2 
TS aa ae yi expl—s xi (p2—e+ 952) 72] 
and this is the appropriate generalization of [77] for two layers. 

We see at once that if the layers are close together (i.e., if their distance of 
separation < \(R/h)'/*) they behave like a singl® layer whose (Am) equals 
(Am),+(Am)>. Having regard to the condition on & {x?/3(e—369)} for impor- 
tant modes sanisionsa above, it is easy to verify that such close lavers make 
p > 1, in agreement with our hypothesis. 

If 

p(Am)> 


(“J OS Evia <P in)” 


the effect of the upper layer becomes negligible (equation [6]). Now the right 
of [7] 
_ (Am) (axst 11 


-—— aren aon 1s \ 3/2 / eae 
= ; x [2.X1 = 20 —(po- 560 , 
(Am), pxettn)" explsy t (or e+ 361) Yi — (p2—€ + 362)" ¥2 | 


and the modulus of this 


(Am)> Bo 
* ay exp{ —2V/2.x(v1V/p1— ¥2V 2). S (€)} 


provided that the layer separation is small compared with the height of the 
lowest laver. Taking y: and y2 approximately equal to unity (mm, m2 are small) 
we see at once that [7] will be satisfied only if* 


[8] —2V/2.x(Vp1— Vp2). F(—€) + log 


The greatest value —x*/?.f(e) can have is odie not greater than 3. It 
is certainly not greater than the airless earth value of 1.606. Thus [7] will 
certainly be satisfied if 


[9] Ah <j RVI tog & 


(Am): 
where Ah is the layer separation, h is the height of the lowest layer, R is the 


*If [7] is satisfied our hypothesis p>> 1 is automatically satisfied. 
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earth’s radius, and d is the wavelength. In deriving this from [8] we have used 
the fact that, for high level inversions, h > R1/8 \?/8, 

We see at once that a first requirement for the lower layer to have the 
dominant effect on the propagation is that (Am)./(A4m), should be fairly 
large (say >e?): the separation must not then exceed the order of R!/® 1/3 h!/?, 
For high level inversions this last quantity is always small compared with h. 
Hence we can say that, for the lower layer to dominate the propagation, 


(i) its (Am) should be considerably greater than the (Am) of the top layer, 
(ii) the layer separation should be small compared with R!/® \!/3 h}/? (i.e. 
the layers are ‘‘close together’). 


If, however, 


[10] 1? octane i 


the upper layer has the predominant effect upon the propagation. Reversing 
the argument given above, we find [10] to be equivalent to 
Am)» 

11 Ah > sRrryinyns log SA22 ab 

This is always satisfied if (Am),/(A4m)2 > e? and also when Ah is large com- 
pared with R!/6 1/3 h!/2, We have seen that this quantity is always less than 
the lower layer height 4, for high level inversions, and, in the case of centimetric 
waves, it is very small compared with 4. Hence we may say that for “widely 


separated” layers (i.e. layers whose separation is not very small compared 
with the height of the lowest layer), it is the upper layer which exerts the 
predominating effect upon the propagation.* 


Note added in proof: \t is of interest to note that Miller (J. Appl. Phys. 22: 
55. 1951) arrives at the same conclusion as that of the present writer’s Part I 
(Can. J. Phys. 33: 241. 1955), using a different technique. 


*We are making the reasonable assumption that (Am); and (Am)2 are of a similar order of 
magnitude. 





NOTES 


PERFORMANCE OF A FAST NEUTRON COLLIMATOR 
By J. T. Sampe, G. C. NEILsON, AND J. B. WARREN 


In the course of an attempt to measure fast neutron polarization by small 
angle scattering in lead (1), a neutron collimator, similar to that mentioned 
by Segel, Schwartz, and Owen (2), was constructed (Fig. 1). 

The experimental arrangement for investigating the degree of collimation 
achieved is shown in Fig. 2. Pulses from a 2 cm. by 2. cm. by 4 cm. deep stilbene 
crystal mounted on a 1P21 photomultiplier were amplified and displayed on a 
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Fic. 1. Horizontal and vertical sections of the collimator. The collimation channel is 
} in. by 1 in., the baffle holes are 3 in. in diameter. 
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Fic. 2. Experimental arrangement for investigating collimation. 
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30 channel kicksorter. At each angle the background from capture gamma 
rays and room-scattered neutrons was determined by inserting a lucite plug 
in the collimator. 

By choosing a suitable pulse height interval between noise pulses and pulses 
due to 2.2 Mev. capture gamma rays from the collimator, the ratio of colli- 
mated neutron counts to background was made 8.7 : 1. The counting rate for 
this interval, less background, is shown in Fig. 3 as a function of angle for 
neutrons of 3.4 Mev. While good collimation was achieved, the number of 
neutrons at angles larger than 3° was too large for the experiment under 
consideration. 
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Fic. 3. Angular distribution of collimated neutrons. 


The authors wish to thank Mr. J. B. Elliott for aid in construction and 
measurement. This work was made possible by financial assistance from Atomic 
Energy of Canada Limited. 
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SOME OBSERVATIONS OF SURFACE PHENOMENA APPEARING ON 
LEAD SINGLE CRYSTALS 


By H. A. Atwater, A. R. LANG, AND B. CHALMERS 


In the course of preparation of some single crystals of lead according to the 
method of Chalmers (2), the occurrence of {111} and {100} crystallographic 
surfaces under certain conditions of freezing was observed. Single crystals of 
lead were grown in a horizontal graphite boat in air, with the oxide film being 
removed occasionally by scraping the liquid surface with a preheated glass 
rod. By using care to introduce a minimum of thermal and mechanical dis- 
turbance, a crystal with very little oxide film could be obtained. The lead used 
was stated to be of 99.999+% purity. Two growth procedures were employed: 
in the first, the boat was suspended so as to reduce mechanical vibration of the 
melt to a minimum; in the second, constant vibration was applied to the boat 
so as to freeze a solid having a rippled surface. 

In the first type of experiment a characteristic step-like structure forms on 
the top surface of the crystal if a {111} plane lies slightly out of coincidence 
with the mean top surface of the melt. The upper surface of the solid at the 
interface tends to grow forward parallel to a {111} plane during advance of the 
interface until further extension is limited by sudden formation of a step in the 
surface. The extending {111} plane tends either to rise above or to penetrate 
below the free liquid surface, depending on its angle with the horizontal. The 
height of the liquid meniscus is appreciably greater than the height of the 
steps which are formed. This is illustrated schematically in Figs. 1(a) and 1(d). 


LIQUID 


(b) 


Fic. 1. Advancing interface with surface steps forming (not to scale). 


In the second type of experiment, a constant vibration of low amplitude 
was applied to the melt. As a consequence of this, a dense pattern of ripples 
approximately 0.2 mm. wide and transverse to the boat axis was frozen into 
the crystal surface. The surface so prepared was found to have highly direc- 
tional reflectivity for visible light, with reflection maxima appreciably sharper 
than those from etch pits in the same metal. The poles of the planes producing 
the reflections were located by means of an optical goniometer of the type 
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Fic. 2. Stereographic projections of two crystals, showing crystal axes obtained by, X-ray 
diffraction (solid points) and axes of surface planes producing optical reflections (open points). 
Arrow at G shows growth direction. 
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Fic. 3. Micrograph of rippled surface (85 X). 
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described by Barrett (1). The poles thus obtained are shown superposed on the 
stereographic projection of the corresponding crystal region in Figs. 2(a) and 
2(6), for two different crystal surfaces. (The axis of the containing boat is 
horizontal in these diagrams, and growth direction was right to left.) 

In each case, optical reflection maxima were found corresponding to re- 
flections from those {111} planes of the crystal which lay most nearly parallel 
to the local top surface. In case 1(a), one of the {100} planes was also thus 
favorably oriented, and produced a strong reflection as well. The angular 
resolution of the reflections was good enough for striations (4) occurring in 
case 1(a) to be readily separated in angle. 

Under magnification, the frozen ripple surfaces appear to be composed of 
segments of plane facets. Fig. 3 shows an example of strong faceting on the 
crystal of 2(a) (magnification 85 X). The conspicuous plane segments in this 
figure correspond to {111} planes on the side of the ripple toward the melt. 
The plane facets on the side of the ripple away from the melt are less extensive 
in area, both for the {100} planes in case 2(a) and for the {111} planes in case 
2(d). 

The occurrence of characteristic surface structures arising in the freezing 
of a metal has also been reported by Elbaum and Chalmers (3), who found that 
the interface exposed by decanting the liquid from freezing lead single crystals 
exhibited a step-like structure when a {111} or {100} crystal plane was within 
approximately 15° of the mear plane of the interface. 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Fine Structure of NO* and NO Emission Spectra in the Schumann Region* 


The recentiy found bands of NO* in the Schumann region (Refs. 1, 5) were photographed 
in the first and fourth order of a 10 ft. vacuum grating spectrograph (dispersion 2.5 and 0.62 
A/mm. respectively). A hollow cathode discharge in fast streaming helium with a small ad- 
mixture of NO gas was used as a light source. The arrangement with a cylindrical LiF lens 
described by Brix and Herzberg (2) acted as a coarse monochromator in the case of the fourth 
order pictures for which the blaze of the grating was highly efficient. The spectra show only 
traces of nitrogen bands in the Schumann region. The partly overlapping NO* and NO (8’) 
bands in the region Ad 2000-1350 A are well developed on Ilford Q2 plates in one-hour ex- 
posures. Second order iron arc lines served as standards for the wavelength measurements. 

NO*.—Table I contains the constants of the NO* molecule resulting from fine structure 
analyses of the following seven bands on the fourth order plates: 1383.0 (1,0); 1413.5 (0.1); 
1428.7 (1,2); 1461.2 (0,2); 1511.5 (0,3); 1564.5 (0,4); 1620.8 (0,5) A. Since the rotational 


TABLE I 
CONSTANTS OF NEW LEVELS OF NO* ann NO (cm.~) 


ZF; A B, a D r,(A) We WeXe 
NOt 
All  73469.6 1.587 0.024 8 X10-§ 1.1926 1608.9 23.3 
i ee 0 2.002 0.0202 6X10-§ 1.0619 2377.1 16.35 
BA; 60365 -2 1.330 0.019 _— 1.303 1216.2 16.4 


structure was fully resolved, the analysis was straightforward. Simple P, Q, and R branches 
could be picked out by inspection. As expected the NO* system revealed itself as a 'II-'2* 
transition in complete analogy with the well-known transitions in the isoelectronic molecules 
N: and CO, and the previously given vibrational analysis was fully confirmed. A small rota- 
tional perturbation occurs at J = 24inthev’ = O level. The internuclear distance in the ground 
state of NO*+ comes out 3% smaller than in the Nz molecule owing to the increased nuclear 
——- of one component of the molecule. 

NO.—Some preliminary results of the rotational structure analyses of the four bands 
1633.5 (1,0); 1739.3 (1,2); 1775.9 (0,2); 1835.1 (0,3) in the 8’ system of the neutral NO molecule 
are also included in Table I. Under low dispersion (Refs. 1,6) the bands show two double- 
headed subbands whose separations agree with the well-known doublet — of the *II, 
ground state. The fourth order plates reveal a rather complex fine structure. Each band consists 
of 12 branches of comparable intensity, which could be followed by calculating the combination 
differences with the help of the tables (Ref. 4) of the ground state rotational levels of both the 
*11, and the *I1,y doublet components. Contrary to the interpretation as *Z-*II given on the 
basis of low dispersion plates, the system is now ascribed to a *A;—-*II, transition with the 74; 
state not far from case b coupling. For low J values, one has rather strong Ra, Qa, Pa (7414-7113) 
and Ru, Que Pu ~~ branches for which AZ = 0, and for higher J a case b) strong 
Re, Qz2, P22, Ru, Qu, Pu branches for which AK = AJ. Using the Hill - Van Vleck formula, 
one obtains A = —2 cm.—!. A small value for the coupling constant A is in accordance with 
Mulliken’s prediction for the 2A state resulting from the... (xo) (wx)* (vx)? electron configu- 
ration proposed earlier for the B’ level (Ref. 1); also, the new assignment of B’ asa A state is in 
agreement with the fact that the B’- BI system (Refs. 1, 3) in the visible shows double double- 
headed bands, and finally that bands corresponding to the B’- A*E transition have never been 


observed. 
More detailed papers on the NOt and NO systems will be published in Helvetica Physica Acta. 


The experimental part of this work was carried out during my stay at the Division of Physics 
of the National Research Council of Canada. I wish to thank Dr. G. Herzberg for his kind 
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hospitality and helpful discussions, and also Dr. A. E. Douglas and Mr. J. Shoosmith for 
valuable help during my stay at the National Research Council in Ottawa. 
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(i) General. Manuscripts, in English or French, should be typewritten, double 
spaced, on paper 8} X11 in. The original and one copy are to be submitted. 
Tables and captions for the figures should be placed at the end of the manuscript. 
Every sheet of the manuscript should be numbered. 

Style, arrangement, spelling, and abbreviations should conform to the usage of 
this journal. Names of all simple compounds, rather than their formulas, should be 
used in the text. Greek letters or unusual signs should be written plainly or explained 
by marginal notes. Superscripts and subscripts must be legible and carefully placed. 

Manuscripts should be carefully checked before they are submitted; authors will 
be charged for changes made in the proof that are considered excessive. 

(ii) Abstract. An abstract of not more than about 200 words, indicating the 
scope of the work and the principal findings, is required, except in Notes. 

(iii) References. References should be listed alphabetically by authors’ 
names, numbered, and typed after the text. The form of the citations should be 
that used in this journal; in references to papers in periodicals, titles should not 
be given and only initial page numbers are required. All citations should be checked 
with the original articles and each one referred to in the text by the key number. 

(iv) Tables. Tables should be numbered in roman numerals and each table 
referred to in the text. Titles should always be given but should be brief; column 
headings should be brief and descriptive matter in the tables confined to a minimum. 
Vertical rules should be used only when they are essential. Numerous small tables 
should be avoided. 


Illustrations 

(i) General. All figures (including each figure of the plates) should be num- 
bered consecutively from 1 up, in arabic numerals, and each figure referred to in the 
text. The author’s name, title of the paper, and figure number should be written in 
the lower left corner of the sheets on which the illustrations appear. Captions should 
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(ii) Line Drawings. Drawings should be carefully made with India ink on 
white drawing paper, blue tracing linen, or co-ordinate paper ruled in blue only; 
any co-ordinate lines that are to appear in the reproduction should be ruled in black 
ink. Paper ruled in green, yellow, or red should not be used unless it is desired to 
have all the co-ordinate lines show. All lines should be of sufficient thickness to 
reproduce well. Decimal points, periods, and stippled dots should be solid black 
circles large enough to be reduced if necessary. Letters and numerals should be 
neatly made, preferably with a stencil (do NOT use typewriting) and be of such 
size that the smallest lettering will be not less than 1 mm. high when reproduced 
in a cut 3 in. wide. 

Many drawings are made too large; originals should not be more than 2 or 3 
times the size of the desired reproduction. In large drawings or groups of drawings 
the ratio of height to width should conform to that of a journal page but the height 
should be adjusted to make allowance for the caption. 

The original drawings and one set of clear copies (e.g. small photographs) 
are to be submitted. 

(iii) Photographs. Prints should be made on glossy paper, with strong con- 
trasts. They should be trimimed so that essential features only are shown and mounted 
carefully, with rubber cement, on white cardboard with no space or only a very small 
space (less than 1 mm.) between them. In mounting, full use of the space available 
should be made (to reduce the number of cuts required) and the ratio of height 
to width should correspond to that of a journal page (4}X7} in.); however, 
allowance must be made for the captions. Photographs or groups of photographs 
should not be more than 2 or 3 times the size of the desired reproduction. 

Photographs are to be submitted in duplicate; if they are to be reproduced 
in groups one set should be mounted, the duplicate set unmounted. 
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Charges for reprints are based on the number of printed pages, which may be 
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requisition form must be ordered officially as soon as the paper has been accepted for 
publication. 





Contents 


The Velocity of Propagation of Electromagnetic Waves—J. T. Hen- 
derson and A. G. Mungall - - - - - - +--+ - = - = 


An Apparatus for Comparison of Thermocouples—T. M. Dauphinee 


An Investigation of the Nuclear Resonance Absorption Spectrum 
of Al?’ in a Single Crystal of Euclase—R. G. Eades- - - - 


Masses of Nuclei in Isobaric Sequences—Kailash Kumar and M. A. 
Preston - - - = -+ 2+ = 2+ 22s 2 = © = © © «= «© 


Long Distance V.H.F. Fields. II. Refractivity Profiles Containing 
‘“‘Sharp Layers’’—F. H. Northover - - - - = - = = = 


Long Distance V.H.F. Fields. III. The Case of Two Elevated Layers 
—Francis H. Northover- - - - - +--+ +--+ e = 


Performance of a Fast Neutron Collimator—J. T. Sample, G. C. 
Neilson, and J.B. Warren - - - = = © © «= © © «= -« 


Some Observations of Surface Phenomena Appearing on Lead Single 
Crystals—H. A. Atwater, A. R. Lang, and B. Chalmers - - 


Letter to the Editor: 


Fine Structure of NO* and NO Emission Spectra in the Schu- 
mann Region—E. Miescher - - - - - = - = = = = 


Printed in Canada by University of Toronto Press 


Page 


265 


275 





